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CALCULATION OF THE LONGITUDINAL AERODYNAMIC 


CHARACTERISTICS OF WING-FLAP CONFIGURATIONS 
WITH EXTERNALLY BLOWN FLAPS 

by Michael R. Mendenhall, Selden B. Spangler, 

Jack N. Nielsen, and Frederick K. Goodwin 
Nielsen Engineering & Research, Inc, 

SUMMARY 

A theoretical investigation has been carried out to extend and 
improve an existing method for predicting the longitucJinal aerodynamic 
characteristics of wing-flap configurations with externally blown flaps 
(EBF) • Two potential flow models were incorporated into the prediction 
method: a wing and flap lifting-surface model and a turbofan engine 

wake model. The wing-flap model uses a vortex-lattice approach to repre- 
sent the wing and flaps. The wing may have an arbitrary planform and 
camber and twist. The flap system may be made up of multiple elements 
of arbitrary size which can have large deflection angles. The jet wake 
model consists of a series of closely spaced vortex rings normal to a 
centerline which may have vertical and lateral curvature to conform to 
the local flow field beneath the wing and flaps. The vortex ring wake 
model may be circular or elliptic in cross section and the wake can be 
specified such that the mass, momentum, and spreading rates are similar 
to the actual turbulent wake. Measured turbulent jet spreading rates 
may be used to specify the jet if analytical spreading rates are not 
available. An iteration procedure is provided to locate the engine 
wake centerline (s) . 

Comparisons of measured and predicted pressure distributions, span- 
load distributions on each lifting surface, and total lift and pitching- 
moment coefficients on swept and unswept EBF configurations are included 
in this report. A wide range of thrust coefficients and flap deflection 
angles is considered at angles of attack up to the onset of stall. 

Results indicate that overall lift and pitching-moment coefficients are 
predicted reasonably well over the entire range. The predicted detailed 
load distributions are qualitatively correct and show the peaked loads 
at the jet impingement points, but the widths and heights of the load 
peaks are not consistently predicted. 



INTRODUCTION 




The short take-off and landing requirement for.STOL aircraft neces- 
sitates a means of achieving very high lift coefficients on aircraft in 
take-off or landing configuration with little sacrifice in cruise 
performance. The externally blown jet-augmented flap provides such a 
means. The jet efflux from engines mounted beneath the wing is allowed 
to impinge directly on the slotted flap system, thus producing a large 
amount of additional lift through engine wake deflection and mutual 
interference effects. 

Several analytical methods for predicting the longitudinal aero- 
dynamic characteristics of EBF configurations have been developed. Some 
of these methods require extensive use of empirically determined infor- 
mation (ref. 1) , and others require little empirical data as input 
(refs. 2 and 3). Reference 3 involves the use of two potential flow 
models; a wing-flap, lifting-surface model and a jet wake model. The 
wing-flap, lifting-surface model is a nonplanar, nonlinear, vortex- 
lattice scheme. The jet wake model consists of a series of vortex ring 
singularities representing a turbulent jet, the boundary of which may 
be specified analytically or with the use of empirical data. These two 
flow models are combined in such a fashion that the tangency boundary 
condition is satisfied on the wing and flap surfaces. Jet wake deflection 
is accommodated through calculation of the flow field beneath the wing 
and estimates of the resulting wake deflection by the user of the method. 

The purpose of the work described herein is to make certain improve- 
ments and extensions to the original prediction method. These fall into 
three categories. In the first, the vortex-lattice model of the wing-flap 
is made more versatile and accurate for modeling multiple flap arrange- 
ments. In the second, the model of the jet wake is extended to include 
elliptical as well as circular cross sections in order to more accurately 
represent the wake characteristics near the impingement area on the 
flaps. In the third, the method of determining the wake centerline 
position is refined and automated within the computer program. 

This report presents the technical approach to the problem, the 
development of the flow models, and comparisons with data. The discussion 
assumes some familiarity with the earlier work of reference 3. A user's 
manual for the resulting computer program is contained in reference 4. 
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SYMBOLS 


a semi-major axis of ellipse 

a^ ,b imaginary and real parts, respectively, of complex roots of 

^ quadratic equation. Appendix B 


initial area of jet wake 


b 


semi-minor axis of ellipse 


c 


chord of area element on wing or flap; wing chord; or 
distance from origin of ellipse and focus, figure 6(b) 


cd u, cn u, 
dn u, sn u, 
tn u 




C 


M- 


d 

E(k) 


f ,f ,f 
u’ v’ w 


Jacobian elliptic functions. Appendix B 


section normal-force coefficient 

axial-force coefficient on an area element on wing, axial 
force/qS; positive as shown in figure 4(a) 

axial-force coefficient on an area element on a flap, axial 
force/qS; positive as shown in figure 4(b) 

drag coefficient, drag/qS; positive aft 

lift coefficient, lift/qS 

pitching-moment coefficient, moment/ qS ; positive in direction 
of angle of attack 

normal-force coefficient on an area element on wing, normal 
force/qS; positive as shown in figure 4(a) 

normal-force coefficient on an area element on a flap, 
normal force/qS; positive as shown in figure 4(b) 

thrust coefficient of a single jet, thrust/ qS ; positive in 
forward direction 

side-force coefficient on an area element on a flap, side 
force/qS; positive as shown in figure 4(b) 

total thrust coefficient for a configuration with multiple 
jets, thrust/ qS 

diameter of circular cross-section jet wake 

elliptic integral of the first kind with argument k 

influence coefficients for a vortex filament, equations (16) 
and (20); positive as shown in figute 3 
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F(k) 


elliptic integral of the second kind with argument k 


'N 


influence function for velocity normal to camber surface, 
equation (3) 




influence coefficients for a wing horseshoe vortex, positive 
as shown in figure 2(a) 



i, j,k 


k 

k' 


influence coefficients for a flap horseshoe vortex, positive 
as shown in figure 2(b) 

unit vectors in x,y,z directions, respectively, see 
figure 6 

integrals defined by equations (A-24) and (B-1) 

argument of elliptic integrals, see equations (32) and 
(B-10) 

1 - k^ 


M 

NF( j) 

NFLAPS 

P 

q 

r 

R 


s 


S 

t 

u,v,w 


U^, v^,w^ 


u . 


, V . , w . 


total mamber of vortex area elements on wing surface 
niomber of vortex area elements on flap j 
total number of flap surfaces 
perimeter of a vortex ring, equation (31) 
free-stream dynamic pressure 

radial distance from a circular jet centerline to a field 
point 

radial distance from a point on a vortex ring to a field 
point, equation (A-15) ; or local radius of a circular jet 

semispan of a horseshoe vortex measured in wing or flap 
chordal plane; or distance measured along jet centerline, 
figure 7 

wing area 

transformed coordinate t), equation (A-16) 

velocities induced at the midpoint of the bound leg of a 
wing horseshoe vortex in the X,Y,Z directions 

velocities induced at the midpoint of the bound leg of a 
flap horseshoe vortex in the X^,Y^,Z^ directions 

velocities in the X,Y,Z directions induced by the jet 
wake 


u 


n- 


V , w 
n’ r 


total induced velocity field at a point on the jet wake 
centerline, equation (34) 


V 


1 



V 


y > z 


y-direction velocity induced at the three-fourth chord 
location of the side edge of a wing area element, see 
equation (24) 

velocities induced at the three-fourth chord location of 
the side edge of a flap area element in the Yf and Zf 
directions, respectively 

vortex ring induced velocities, equations (A-25) , (A-26) , 

and (A-27) 


V 


x,y,z 




X,Y,Z 


free-stream velocity 

initial jet wake velocity, equation (29) 

coordinate system fixed at center of ellipse, figure 6; or 
coordinate system associated with a lifting surface 
containing a horseshoe vortex 

jet coordinate system fixed at the center of the jet inlet, 
figure 7 

wing coordinate system, see figure 2(a), with origin at 
wing root-chord leading edge 


Xf,Yf,Zf flap coordinate system, see figure 2(b), with origin at 

inboard side-edge leading edge of flap 


a 


wing root chord angle of attack relative to free-stream 
direction 




angle between tangent to mean camber surface of the wing 
and the X direction 


Y 

r 



e 


0 

p 

0 


jet wake vortex cylinder strength 

vortex ring strength, or horseshoe vortex strength 

angle between tangent to mean camber surface of a flap and 
the Xf direction 

streamwise flap deflection angle measured in a plane 
parallel to the X-Z plane, positive downwards 

pressure difference across a lifting surface 

flow angle of jet centerline measured in x.,z. plane, 
equation (35) ^ ^ 

average flow angle of jet centerline, equation (37) 

elliptic coordinate. Appendix A 

density 

flow angle of jet centerline measured in x. ,y- plane, 
equation (36) 
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a 


0 


0f 


f 

j 

left 

o 

right 

S 

W or WF 
a 

I, II, III, 
IV 

1/2 


average lateral flow angle of centerline, equation (40) 
constant defined in equation (B-39) 

wing dihedral angle, positive as shown in figure 2(a) 

flap dihedral angle, positive as shown in figure 2(b) 

sweep angle of the bound leg of a wing horseshoe vortex 
measured in the wing chordal plane, positive swept back 

sweep angle of the bound leg of a flap horseshoe vortex 
measured in the flow chordal plane, positive swept back 

Subscripts 

flap 
jet wake 

left leg of horseshoe vortex 

initial conditions 

right leg of horseshoe vortex 

slat 

wing-flap 
free stream 

denote four quadrants of elliptic vortex ring, figure 6(a) 

denotes position at which jet axial velocity is half the 
maximvun velocity on the centerline 

Superscripts 

vector 

primes indicate local coordinate system 


6 



ANALYSIS 


The work described in this report is basically an extension of that 
reported in reference 3 . Certain improvements were made in both the 
lifting-surface model and the jet model, and both models were then 
combined to form an iterative approach to the calculation of the aero- 
dynamic characteristics of EBF configurations. In this section, each 
flow model and its modifications are discussed individually, and then 
the total EBF interference procedure is described. 

Wing-Flap vortex -Lattice Model 

General description .- The vortex-lattice lifting-surface model used 
in the present investigation is an extension of that used in references 3 
and 5. In reference 3, the configuration treated consisted of a wing 
with a single trailing-edge flap with a uniform deflection angle. 

Multiple flaps were handled by specifying a camber distribution over 
some mean flap surface. In reference 5, multiple spanwise flap segments 
were considered, but separate loading distributions on individual flap 
elements were not available. The vortex-lattice model in both refer- 
ences 3 and 5 allowed the trailing vortex legs to extend aft to infinity 
in the plane of the originating surface; that is, the wing chordal plane 
or the flap chordal plane. 

Typical externally blown flap (EBP) and upper-surf ace-blowing (USB) 
configurations have multiple spanwise flap segments, each of which may 
be made up of multiple flap elements. in the present investigation, the 
method and associated computer program have been modified so that each 
flap segment is defined as a separate lifting surface for the purpose of 
computing forces and moments. The trailing vorticity is deflected so 
that it lies on each surface aft of its originating surface and then 
extends to infinity in the plane of the last surface as illustrated in 
figure 1. 

The wing and flap configuration parameters included in the method 
are listed below. 

Wing Panels 

Leading-edge shape: May have up to 30 breaks in sweep. 

Trailing-edge shape: Same as for leading edge. 

Taper: Determined from leading-edge and trailing-edge shapes. 
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Tip chord: parallel to root chord. 

Dihedral: Arbitrary but constant over the semispan. 

Mean camber surface: May have both twist and camber. 

Thickness: Neglected. 

Flap Panels 

Number: Up to 10 individual flap segments; maximiam of 3 at any 

spanwise region. 

Location: First flap in any spanwise region is at or near wing 

trailing edge; gaps are permitted between wing and first flap and 
succeeding flaps. 

Leading-edge shape: Straight line which may be swept. 

Trailing-edge shape: Same as for leading edge. 

Taper: Linear. 

Root chord: In plane parallel to vertical plane containing wing 

root chord . 

Tip chord: Parallel to root chord. 

Span: Full or partial span. 

Deflection: Each flap may have a different deflection angle. 

Mean camber surface: Each flap may have both twist and camber. 

Thickness: Neglected. 

The method is restricted to angle of attack; that is, yaw effects 
are not included. Compressibility effects are not included since large 
flap deflections are only used at low speeds. Angle of attack, flap 
deflection angle, dihedral angle, and twist and camber angle are accounted 
for in the boundary condition by trigonometric functions since some of 
these angles may be large. 

The vortex-lattice arrangement and the coordinate system for a 
swept wing with trailing-edge flaps deflected are shown in figure 1. The 
wing and flaps are divided into trapezoidal area elements. The flap side 
edges must lie along chordwise edges of the area elements, and the span- 
wise distribution of area elements on the flaps must be the same as on 
the wing section ahead of the flap. 
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A horseshoe vortex is placed in each area element such that the 
bound leg lies along the quarter chord of the, area element and the 
trailing legs lie along the side edges of the element in the plane of 
the element. The trailing legs trail back in the plane of the area 
element until the leading edge of the next surface is encountered. They 
are then deflected so that they lie in the plane of that surface. This 
deflection is continued through the last surface. The trailing legs 
are then extended aft to infinity in the plane of the last surface. 

The flow tangency boundary condition is applied at the midspan of 
the three-quarter chord line of each area element. The wing and flap 
boundary conditions will now be determined. 

Wing boundary condition .- The wing boundary condition is that the 
velocity normal to the mean camber surface at a wing control point be 
zero. This can be written 


M 

NFLAPS 

NF(j) . 


Z Vv.n 

n=l 

^ I 

j=l 

) + V. 1=0 

Z_ r,p,v,n N,V 

n=l ^ ' 

(1) 


V = 1, 2, . . .M 

The first term is the sum of the normal velocities induced at the 
control point by the M wing vortices. The second term is the sum of 
the normal velocities induced by all of the flap vortices. There are 
NFLAPS flaps and NF(j) vortices on the flap. The component of the 

free-stream velocity normal to the camber surface is the third term, 

The last term contains externally induced velocities such as those 
due to an engine wake. The derivation of the four terms will now be 
presented . 

Figure 2 (a) shows coordinate systems associated with a wing inclined 
at dihedral angle, 0. The X,Y,Z system is the basic wing system shown 
in figure 1. The quantities F,^,F^,F,^ are influence function^, dimen- 
sionless velocities per unit vortex strength, with positive directions 
as shown. The X',Y',Z' system is obtained by rotating the X,Y,Z 
system through the wing dihedral angle, 0. The quantities F^,F^,F^ 
are influence functions in the primed coordinate system and are 
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( 2 ) 


F ' = F 
u u 




F' = F cos 0 + F sin 0 / 

■V -V V ' 


F' = -F sin 0 + F cos 0 
w V w 


The 


Consider the following sketch which is drawn in the Y' = 0 plane, 
angle is the local "angle of attack" due to twist and camber. 



The influence function normal to the surface is 


N 


F^ cos sin a_g 


(3) 


If Vp ^ ^ is the normal velocity induced by the n wing vortex at 
^ w’ ’ 


the control point 


V 


^n 


V 


47tV n 


(4) 


v,n 


Then, using equations (2) and (3), the first term in equation (1) is 


M 


M 


^ V' V, = V ( F cos 0 cos a. - F sin 0 cos a 

V ^ r^,v,n 47TV ^ w^^^ '"v v^_^ '"v 


n=l 


n=l 


v,n 


+ F sin a„ \ 

^,n K) 


(5) 


In order to evaluate the second term in equation (1) , the wing 
boundary condition, the flap coordinate system must be defined. The wing 
and flap coordinate systems are shown in figure 2 (b) . The angle 6 is 
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the flap deflection angle and is measured in a strecimwise plane parallel 
to the X,Z plane. The flap coordinate system is obtained by rotating 
the wing system about the Y axis by the angle 6^^ and then translating 
to rhe flap inboard-side leading edge. The influence functions 
and F are known in the flap coordinate system. They can be resolved 

Wf 

into the wing coordinate system by 




F = F cos 6 + F sin 5 „ 

u u^ xz w^ xz 


F = F 


F = -F sin 6 + F cos 6 

w u^ xz w^ xz 


( 6 ) 


Using the above expressions and equations (2) and (3) , the influence 
function normal to the local camber surface is 

F„ = F (cos 0 cos a^cos 6 + sin a.sin 6 ) - F sin 0 cos a„ 

w^' ^ £ xz £ xz' ^ 

+ F (sin a. cos 6 - cos 0 cos a.sin 6 ) (7) 

u^ i xz ^ i xz' 

The normal velocity induced by the n^^ flap vortex at the wing 

control point is 


V 


r^TjV.n r 


V 


4ttV N^,v,n 


( 8 ) 


The second term in equation (1) is evaluated by siamming over all flap 
vortices 


NFLAPS 

NF(j) 

NFLAPS 

NF(j) 

v" I 

/ ■ r o > V j n 

- Z 

Z 

n=l 

n=l ^ 

j = l 

n=l 



+ sin 

sin 6 


+ F /sina„ cos 6 „ - cos 0 cos a. sin 6 ) 

u^,v,n{ £^ xz^ £^ xZj^J 


(9) 
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The third term in equation (1) is the component of the free-stream 
velocity normal to the wing camber surface at the v ” control point. 

X 



From the above sketch it is seen that 

F = -V cos a 
u 

V / 

F = -V sin a 
w 

Equations (2)., (3) , and (10) give 


( 10 ) 


^ = -sin a cos 0 cos a. - cos a sin a. (11) 

The fourth term in equation (4) is the velocity normal to the wing 
surface produced by other sources. in the present application, these 
sources are the engine wakes. Assume these velocities are known in the 
wing coordinate system and are positive in the positive X,Y,Z directions 
Then 


F = u . 

U X, V 


F = V. > 
V 1, V ^ 


F = W . 

W X, V 


th 


at the V control point, and with equations (2) and (3) , 
V. 


w . 

V V 


V . u . 

cos 0 cos a. - — Y* ^ sin 0 cos a. + ■ sin a 


^ V 


V 


^ ■ V 


(12) 


(13) 
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Substituting equations (5) , (9), (11), and (13) into equation (1) 

and putting the known quantities on the right-hand side of the equal 
sign, the final form of the wing boundary condition at the v control 
point on the wing is 


M 


V* F cos 0 cosOfl -F sin 0 cos a „ +P sina, 

L. 4ttV I. W V u^^^ 

n=l ’ 

NFLAPS NF(j) p 

+ > -r-^ F fcos 0 cosa« cos 6 +sina. sin 

^ Z-. 4ttV ''w,,v,n', xz^ xz^; 

j=l n=l 


- F 


sin0 cos a, +F „ ^(sina. cos 6 - cos 0 cos a. sin 6 \ 

:,v,n Up,v,nV xz^ v 


W . V . / ^ • \ 

= ^sina- cos 0^cos + -^ sin 0^cos + (cos ot - ) sin a 


V 

(14) 


Flap boundary conditions .- The flap boundary condition can be 
derived in the same manner as the wing boundary condition. The 
resulting boundary condition at the control point on the i 

flap is 'i 


M 


Z 4^ (cos« cos a cos 6^ -sin 6, sin 6^ ) 

n=l >- v,n V V 


- F 


sin0^ cos 6 - +F /sin 6. cos 6 ^, + cos 0 cos 6 sin 6 ^ 

,n S K '^Vyry \ ^^i 

NFLAPS NF(j) p , 

Z Z 5S; jv.v.n ) 

j=l n=l \ X V V X J 


- sin 6 


„ sin/s -6 \ -F sin0£ cos 6 . 

iv . xz. °xZjjJ Vp,v,n 


+ F si 

Up,v,n|_ 


in S^^cos (6^^ _ - 6^^ _ ) t cos 0^^cos 6^^sin(6^^ _ - 6^^ _ )JJ 

(continued on next page) 
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sin 


(“ * 


cos 0^ cos + cos 




sin 6 



W . V . 

- (cos 0£ COS 6j, cos 6^^ - s±n6. sin 5 ^os 6 - 

'' \ V '*'v i i' V 


u . 

cos 6^2 + cos 0 ^ cos 6 ^ sin 6 ^^ j 


(15) 


In this equation, 0^ is the dihedral angle of the flap containing the 

control point and dp is the local "angle of attack" due to twist 

V 

and camber at this point. The angle °xz • deflection angle of 

the flap containing the control point and . 1® deflection angle 

of the flap in which the influencing vortex originates. 

Influence functions.- The functions F .F .F and Fn .F,, .F,., are 

u’ v’ w ^f’ '"'f 

influence functions relating the velocity components induced at some point 
by a horseshoe vortex to the circulation strength. In reference 3 a 
horseshoe vortex was considered to have three parts; a bound leg and 
two trailing legs extending to infinity in the plane of the surface 
containing the bound leg. Figure 1 shows the horseshoe vortex arrange- 
ment used in the present method. Except for a vortex originating in the 
last surface at a spanwise location, the trailing legs are composed of a 
series of finite length segments followed by two semi-infinite length 
segments. Thus, the equations given in reference 3 for the influence 
functions are not applicable, in general, to the present method. 

General expressions are required for the influence functions for a 
finite length vortex filament and a semi-infinite filament. Such 
expressions are available in reference 6. For completeness, they will 
be repeated here. 

The coordinate system for a finite length vortex segment is shown 
in figure 3(a). The x,y,z coordinate system is that associated with 
the surface in which the horseshoe vortex containing this segment 
originated. In that system there is a vortex filament extending from 
X ,Y ,z to X ,y ,z with vorticity directed as shown. The influence 
functions at point x,y,z are, from the Biot-Savart law, 
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( 16 ) 


( 17 ) 


( 18 ) 


( 19 ) 
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similarly, for the semi-infinite length vortex filament shown in 
figure 3 (b) , 


f 


u 


■P 

“v 



+ f2 


+ f2 



(1 

(1 


COS 0 ' ) *1 


cos 0') ) 


f = 
w 


-f 


>/d2 + 


+ f2 


(^) (1 - cos 0') 


J 


(20) 


where 

d = Ay(z - z^) - A^{y - y^) 

e = A(x-x) - A(z-z) > 
z 1 X 1 ^ 

f = A (y - y ) - A (x - X ) 

X 'J -I 1 ' y ' 1 ' ^ 


(21) 


The positive sense of vorticity is shown in figure 3 (b) . The quantities 
AxjAyjA^, in equations (21) are the direction cosines of the vortex 
filament in a direction from infinity. Also, 


cos 0 ' = 


r ' = 


r -^x^^ ~ + Ay(y - yi) + A^(z - Zi) 1 

L ^ (X - x^)^ + (y - y^)^ + (z - Zi)^ J 

R' ^1 - COS^0' 


^ (x - X^)^ + (y - y^)^ + (z - Z^)' 


J 


( 22 ) 


In order to evaluate the influence functions appearing in equa- 
tions (14) and (15) , equations (16) through (19) are used to evaluate 
the influence of the bound leg and all finite length trailing leg 
segments of both the right and left trailing legs. Equations (20) 
through (22) are used to evaluate the influence of the two semi-infinite 
trailing legs. The influence functions in equations (14) and (15) 
include the influence of the vortex on the left wing panel or left-side 
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flap and that due to the corresponding vortex on the right wing panel 
or right flap at the left-side control point. 

Aerodynamic forces and moments .- Once the strengths of the vortices 
have been determined by solving the set of simultaneous equations given 
by equations (14) and (15) , the aerodynamic forces acting on each area 
element on the wing or flap can be calculated. The Kutta-Joukowski law 
for aerodynamic forces on a vortex filament is used. The method of 
calculating the forces is the same as that described in reference 3 
except that certain terms omitted there are now included. Following the 
procedure of reference 3 the axial-force coefficient produced by one of 
the area elements on the left wing panel is 

(2s) (cos 0 sin a - ^ cos 0 + ^ sin 0) (23) 

Similarly, the normal-force coefficient is 


an 

SV 




left 


right 


(24) 


Positive directions of the forces are shown in figure 4(a) . In the above 
equations, 2s is the width of the panel in the wing chordal plane. The 
velocities u, v, and w, positive in the positive X,Y,Z directions, are 
those induced at the bound-leg midpoint by the entire wing-flap vortex 
lattice plus velocities due to the engine wake, if present. The angle 
0 is the sweep angle of the bound leg, positive if swept back. 

The second term in equation ( 24) is the net force on all of the 
vortex trailing legs lying on the two side edges of the panel. The 
vortices are only those lying in the chordwise strip containing the 
panel. In equation (24), EF is the sum of the strengths of these 
vortices, c is the length of the panel side edge, and v^ is the 
vortex lattice and engine wake induced velocity at the three-quarter 
chord location of the side edge. The siobscripts left and right refer 
to the left side edge and right side edge. 

Similar expressions can be written for the forces acting on a panel 
on a flap. Positive directions of the forces are shown in figure 4(b) . 
The axial-force coefficient is. 
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9r T w- u. 1 

(2s) l^cos 0^cos (a + 6^^;) ~ If T J 


the side-force coefficient is, 

C = ^ 

Yf SV 


r "f 

(2s) cos (a + 6^^)sin (j)^ - -^ sin 0^ - sin (a + 6^^) tan ip^ 




ight J 


(26) 


and the normal-force coefficient is, 


N. 


r V 

= 1^ (2s) tan + cos (a + '5^^)cos 0^ - 


^f . 

— cos 0^ 


SI 

SV 




right -■ 


(27) 


The velocities u^^,v^,w^ are positive in the positive X^,Y^,Z^ 
directions. 

The elemental panel forces given by equations (23) through (27) 
can be used to calculate spanwise load distributions by summing over 
the area elements in a spanwise row and total forces and moments can be 
determined by then summing spanwise. Moments can be determined by using 
the panel forces with the appropriate moment arm. 


Vortex Ring Jet Wake Model 

A potential flow model of the wake from a turbofan engine is 
presented in reference 3. This model has a vorticity distribution placed 
on the specified boundary of an expanding, circular cross-section jet. 
This is modeled by a series of vortex rings coaxial with the prescribed 
centerline of the jet. Each ring represents a finite increment of length 
of the cylinder, and the vortex strength of each individual ring is equal 
to the net vorticity on the incremental length of cylinder which it 
replaces. A sketch of this circular ring model reproduced from refer- 
ence 3, is shown in figure 5 along with some typical velocity profiles 
at various distances downstream of the beginning of the jet. This 
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model represents the mass, momentum, and entrainment characteristics of a 
turbulent, coflowing jet. While the velocity profile within the jet is 
approximated by a uniform profile, the induced velocities outside the 
jet are probably accurate, because they are related to the entrainment 
induced flow. 

Flow surveys taken just aft of the trailing edge of an EBF model 
(ref. 7) indicate that an initially circular jet becomes distorted after 
interaction with the flap system. The distortion is such that the wake 
tends to flatten and become more elliptical than circular in cross 
section. Observations of free, nonaxisymmetric cross-section jets 
(refs. 8 and 9) show that these jets tend to become axisymmetric as they 
progress downstream. Thus, the most probable picture of an engine wake 
is one which has the wake leave the engine as a circular jet, become 
more elliptical as it approaches the flap, flatten as it passes under 
the flap system, and return gradually towards a circular section in the 
far field. In the case of large flap deflection angles, it is possible 
that the wake aft of the flap may roll up into a pair of vortices, as 
has been observed with circular jets issuing normal to a crossflow, but 
there are few observations on which to base any conclusion for EBF 
configurations. 

On the basis of the EBF jet observations noted above, it was felt 
that the method of reference 3 should be augmented with an elliptical 
jet model which can be combined with the circular jet to produce a wake 
which can change continuously from circular to elliptical, with axis 
ratio variation to be arbitrary and selected on the best available 
information. This section, together with Appendices A and B, describe 
the development and use of such a model. 

Velocity induced by an elliptic vortex ring . - Analogous to the 
circular ring, the first requirement of the elliptic jet model is a 
solution for the velocity components induced at an arbitrary field 
point by an elliptic vortex ring. The coordinate system and notation 
are shown in figure 6. 

Unlike the circular ring, no solution was available for the induced 
velocity of the elliptical ring. Consequently, the solution was developed 
and is discussed in detail in Appendix A. Through the use of confocal 
elliptical coordinates, an analytical solution can be obtained. The 
basic solution for the velocity components, for one quadrant in a ring 
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coordinate system, are given as equations (A-23) through (A-25) . These 
equations are used, with appropriate sign changes, to sum the contribu- 
tions of all four quadrants of the ring. These velocity components can 
then be related to components in the jet and wing coordinate system 
through the nomenclature of figure 7 and equations (A-31) and (A-32) . 

The solution is obtained in terms of certain integrals, J^, which are 
developed in Appendix B. 

Elliptic jet model .- The development of the elliptic cross-section 
potential flow jet model proceeds in similar fashion to the development 
of the circular jet model in reference 3. The assumption is made here 
that the jet begins with a circular cross section and becomes elliptical 
at any rate desired. 

Consider for the moment that the jet vortex model is a constant 
radius, semi-infinite circular vortex cylinder with strength y/V. The 
strength of this vortex cylinder is specified by the relation 


I = 

V 



(28) 


where Vj is a uniform velocity across the engine exhaust exit. This 
jet velocity is given approximately as 


V 


1 

2 


1 + 




1 + 2C 




T A . p . 
3^3 


(29) 


where is the thrust coefficient for a single jet and A^ is the 

exhaust exit area for a single jet. 

Starting with a vortex cylinder with constant radius and 

strength y/V, the actual jet wake is modeled by distorting the cylinder 
to fit a prescribed centerline path and outer boundary. The continuous 
vorticity distribution on the new distorted cylinder is replaced by a 
closely spaced series of vortex rings coaxial with the jet centerline. 
Each ring represents a small increment of length of the jet, and each 
ring has a shape which corresponds to the specified jet boundary. The 
strength of the ring vorticity is equal to the net vorticity on the 
incremental length of the original cylinder. Thus, the semi-infinite 
length of the original cylinder is replaced with a finite length of 
vortex rings. As described in reference 3, the semi-infinite length 
cylinder can be approximated to any required accuracy in this manner. 
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The strength of any vortex ring along the jet is 


r 

V 


p 

^ AS 
VP 


( 30 ) 


where P is the perimeter of the ring in question and As is the 

length of the wake represented by the ring. For an ellipse, the perimeter 

is 

P = 4a E(k) (31) 


where E (k) is a complete elliptic integral of the first kind, and 


k^ = 



(32) 


The quantities a and b are the semi-major and semi-minor axes of the 
ellipse. An approximate relation for the perimeter is 

P ^ ^ (33) 

which is within 5 percent of equation (31) for a/b < 3. 

From equation (30) , the strength of the vortex ring decreases as 
its perimeter increases such that the product of r and the ring perime- 
ter is constant for equally spaced rings (As = constant) . This is 
equivalent to keeping the total vorticity on the cylinder constant. The 
induced velocity at any point in the field is calculated by summing the 
contributions from all the rings making up the model. Typical velocity 
profiles inside various elliptic jets are discussed in the results portion 
of this report. 

Jet description .- As is the case in reference 3 for circular jets, 
the current model for elliptic cross-sectional jets requires three items 
to completely determine the analytical description of the jet. The first 
is the initial vortex ring strength from equation (28) , which is related 
to the thrust and momentum in the jet. The remaining two items to be 
specified are the position and path of the centerline of the jet and the 
boundary of the jet. 

The computer program described in reference 4 allows the centerline 
of the jet to move under the influences of the free-stream velocity, the 
wing and flap loading induced flow field, and the jet induced flow field. 
The objective in permitting the centerline to move is to be able to 
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locate the centerline along the streamline of the wing-flap- jet flow 
which leaves the center of the engine exhaust. This is an approximation 
to placing the singularities representing the jet along streamlines, 
which cannot be handled because of the limitations on cross-section 
shape (circle and ellipse) . Thus, the jet centferline is free to move 
so that an iteration can be performed until convergence is attained 
between the centerline and streamline positions. Since each jet is free 
to move according to the local flow field induced along its centerline, 
the user need not specify the jet centerline in great detail prior to a 
calculation. In practical usage, the initial centerline may be assvimed 
to be straight. The motion of the centerline is calculated in the 
following manner. 

The path of the jet centerline is specified by a number of points 
on the centerline. At a point n, the induced flow field has the three 
components 
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are jet induced velocities 
are wing-flap induced velocities 
are free-stream components 


The resultant flow angle in the X-Z plane is 

% - 


and in the X-Y plane is 


(35) 


( 36 ) 
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Thus, the vertical flow angles at two consecutive points along the 
centerline appear as in the sketch below. 



X 


Assuming that the average flow angle between points n-i and n is 
given by 


0 = 


n 


= (0 + 0 ) 

2 n-i n ' 


(37) 


a new vertical coordinate at point 
is given by 


z 

n 


n-i 


n based on the average flow angle 

- X ^tan 0 (38) 

n-i ^ n 


Similarly, a new lateral coordinate at point n is given by 


^n = 


n-i 
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- X ") tan 
n-i / 
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where 


a = 
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(40) 


Starting at the engine exit with fixed coordinates and slope, the entire 
jet centerline is repositioned according to equations (38) and (39) . 

A limit is placed on the vertical deflection angle 0^^ in accordance 
with the following considerations. It is well known experimentally that 
the flaps of an EBF configuration do not achieve 100-percent jet turning 
efficiency but a somewhat lower value which decreases with increasing 
flap angle. It was found in the calculations that with large flap 
deflections (6_„ > 40°) , the predicted jet centerline slope would some- 
times equal or exceed the aft flap deflection angle, with the result 
that the predicted interference effects would be low because the jet was 
passing too far below the flap. Under those circumstances, if a limit 
on 0^ of the product of the measured jet turning efficiency and the 
flap deflection angle were imposed, the jet would pass closer to the 
flap, the interference effects would be larger, and the predicted and 
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measured force results were in good agreement. 'This was felt to be a 
reasonable limitation in view of the experimental turning efficiency- 
results, and was incorporated in the computer program. If the computed 
value of 0^ from equation (37) exceeds a specified limit, 9 ^ is 
set equal to the limiting value. The remainder of the centerline calcu- 
lation and iteration procedure is carried out as before. 

To complete the description of the jet wake, the boundary of the 
jet must be specified at all points along the jet. In reference 3, an 
axisymmetric turbulent jet in a coflowing stream was defined using an 
analytical method presented in reference 10. The radius distribution of 
an axisymmetric jet is shown in figure 8 for various jet velocity ratios. 
Unfortunately, a similar series of curves is not available for noncir- 
cular cross-section jets; therefore, empirical evidence must be used 
where possible. 

Elliptic and rectangular shaped jets exhausting into still air are 
described in references 8 and 9, and measured half-width boundaries are 
shown in figure 9. These data are not sufficient to describe the 
expansion of typical elliptic cross-section wakes occurring on EBF 
configurations, but they are shown here to illustrate several general 
points. The first is that elliptic and rectangular jets with similar 
axis ratios exhibit the same kind of expansion characteristics. It also 
appears that the minor axis expands at nearly the same rate as a circular 
jet. Figure 9 also illustrates that nonaxisyrametric jets exhausting into 
still air tend to become axisymmetric at large distances downstream. 

Measured jet wake velocity profiles at a number of spanwise stations 
approximately one flap chord aft of the trailing edge of the last flap 
of a four-engine EBF configuration are presented in reference 7. The 
boundaries of three jets corresponding to total thrust coefficients of 
1, 2, and 4 were deduced from these data. The thickness distributions 
of the jet are shown in figure 10 as solid lines. The spanwise 
extent of the jets is not well defined by the measurements; therefore, 
the extrapolated inboard and outboard portions of each jet are shown 
dashed. To illustrate the growth and spanwise spreading of each jet, 
the original jet cross-sectional areas are shown at their appropriate 
semispan stations as dashed lines. The relative vertical positions of 
the original jets and the measured profiles are not comparable as the 
jets aft of the flaps have been deflected downward. 
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The results of figure 10 show that the boundary of the jets tends 
to form an elliptical shape after interaction with the flap system. The 
jets also seem to move inboard toward the fuselage instead of outboard 
toward the wing tip. This motion seems to become more pronounced as 
the thrust coefficient increases. 

The data of reference 7 can also be used to examine the expansion 
ratio of the jets. For this purpose, areas of the jets were obtained by 
graphical integration from data like that of figure lO and divided by 
the initial areas of the two jets. The results for the station just 
downstrecim of the flap trailing edge (about 9 initial jet radii down- 
stream of the jet exit plane) are shown in figure 11(a) as a function of 

C and jet velocity ratio (from eq. (29)). The shaded area represents 
M' 

the scatter and the solid line is the least-squares fit to the data. 

If we assume the jet area shown in figure 10 is redistributed into 
a circular cross section for each individual jet, an equivalent radius 
ratio, R/R^, can be calculated for each thrust coefficient as simply 
V A/Aq. This result is shown in figure 11(b) where it is compared with 
the radius ratio of an axisymmetric jet from figure 8. Figure 11(b) 
indicates that an aixisymmetric jet alone expands at a faster rate than 
an initially axisymmetric jet under the influence of a wing and flap. 

This kind of information is important in selecting a spreading rate for 
a typical eBF jet. 

EBF Interference Calculation Procedure 

Calculation of the aerodynamic loading of a wing-flap configuration 
under the influence of the jet wake of a turbofan engine is done with the 
combination of the two potential flow models described in the preceding 
sections. The two flow models are combined by superposition. The jet 
model induces a velocity field on the wing and flap which produces an 
interference loading on the lifting surfaces. The wing and flap loadings 
induce a velocity field in the vicinity of the jet and tend to deflect 
the jet away from these surfaces. It is assiamed that the engine thrust 
is unaffected by the presence of the wing-flap. Because of the mutual 
interaction between the jet and lifting surfaces, either a simultaneous 
solution involving both flow models or an iterative solution is required. 
The iterative approach was chosen in reference 3 and that approach is 
continued, with some refinements, in this report. The solution is carried 
out in the following manner. 
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Before any calculations are made, the jet centerline is positioned 
with respect to the wing and flap. The initial location of the center- 
line can be based on some a priori knowledge of the flow field beneath 
the wing and flap system or it can be located in -a strictly arbitary 
fashion. For example, it is quite acceptable to choose the initial jet 
centerline to be a straight line aft from the engine exhaust. 

The (expanding) jet boundary and cross-section shape distribution 
should be chosen according to whatever procedure seems most appropriate, 
and this distribution, which will be unchanged from this point on, is 
placed on the centerline to define an initial jet wake. The jet induced 
velocity field is computed at selected control points on the lifting 
surfaces and the circulation distribution on the wing and flaps is 
obtained such that the tangency boundary condition is satisfied at each 
control point. 

At this point in the solution, the boundary condition on the wing 
and flap surfaces is satisfied but the jet position has not been influenced 
by the presence of the wing and flap. The wing-flap influence consists 
of modifying the jet location to cause the jet centerline to lie along a 
streamline of the combined jet-wing-flap flow. The initial jet center- 
line is then adjusted by computing the total flow field at a number of 
points on the centerline and moving the centerline to a new position such 
that it lies along the computed flow direction at the specified points. 

This completes the first iteration of the solution. In this situation, 
the wing and flap loading is not compatible with the jet flow field 
corresponding to the new position of the jet. Thus, a second iteration 
is needed. 

The flow field corresponding to the adjusted jet position is computed, 
and a new wing-flap loading distribution is obtained. The jet centerline 
is again moved to lie along the new flow directions. This procedure is 
continued until either the centerline position or the total wing and flap 
loading converges to within a desired tolerance. With a converged 
solution, the total flow is tangent to the wing and flap surfaces and 
the jet centerline lies along a streamline of the flow. 
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RESULTS 


The methods of analysis described in the previous section have been 
evaluated in the following manner. Each flow model (the lifting surface 
and the jet wake) is examined separately to determine its individual 
characteristics. Where possible, comparisons with data are made to 
assist in the evaluation. After the individual flow models are examined, 
they are combined to perform several calculations on EBF configurations. 
Convergence characteristics of the jet wake centerline iteration procedure 
are first examined. Comparisons with experimental data are then made and 
used to evaluate the total EBF prediction method. The data comparisons 
and results are discussed below. 

Wing-Flap Model 

The first configuration used to check out the deflected wake vortex- 
lattice method is an aspect ratio 6 rectangular wing with a single slotted 
flap (ref. 11) . The Fowler-type flap has a chord that is 30 percent of 
the wing-alone chord and is deflected 40°. The vortex-lattice arrangement 
used for this configuration is 8 chordwise by 13 spanwise on the wing 
semispan and 4 by 13 on the flap. 

The predicted results for this wing (and following configurations) 
can be compared with similar results from reference 3 to illustrate two 
features of the present method: the deflected horseshoe vortex trailing 
legs and the gap between the wing and flap. The lattice arrangement is 
somewhat different between the two calculations, but the work of Appendix A 
of reference 3 indicates that the lift and moment should be converged 
with both lattice arrangements. 

Comparisons of measured and predicted lift and pitching-moment 
coefficients are shown in figure 12(a). With zero flap deflection, the 
deflected wake method reduces to the undeflected wake method and the 
results are identical; therefore, only one predicted curve is shown. 

The agreement between measured and predicted lift curves is excellent 
up to the point where the lift curve starts to bend due to separation 
or stall. The predicted pitching -moment curve for this configuration 
is also in reasonably good agreement with the experimental results. 

The effect of the lift on the fuselage is not considered in these 
comparisons. An estimate of the fuselage lift and pitching-moment 
coefficients given by slender-body theory results in values of 0.01 and 
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0.02, respectively, at 10° angle of attack. Inclusion of the fuselage 
pitching moment would tend to improve the agreement on moment curve 
slope. 

The effect of deflecting the trailing legs of the horseshoe 
vortices is illustrated in the upper portion of figure 12(a). In 
general, it is observed that when the trailing vorticity is allowed 
to move aft in the plane of the surface from which it originates, the 
predicted lift coefficient is larger than when the trailing vorticity 
is allowed to bend around the flap surfaces and finally move aft in the 
plane of the last flap. Numerically, the lift coefficients corresponding 
to the undeflected wake are approximately 10 percent higher than those 
from the deflected wake calculations and the lift curve slope is slightly 
higher. The pitching-moment coefficients corresponding to the deflected 
wake are in much better agreement with experiment, both in magnitude and 
slope, than the results obtained considering the wake to be undeflected. 

The model used in reference 11 has a slotted flap, but the size of 
the gap is not specified. The effect of gap size is shown in figure 12(a), 
where predicted results for gaps of 0.02 c and 0.03 c are shown for both 
deflected and undeflected wakes. The results show that increasing the 
gap size decreases the total lift and increases the nose-down pitching 
moment. A closed gap was considered for the previous undeflected wake 
calculations and these results are consistent with those just discussed. 

The effect of the gap size is not insignificant; and when using the 
wing-flap vortex-lattice method described herein, one should take care 
in modeling the gaps between lifting surfaces correctly. 

Measured and predicted chordwise pressure distributions on the wing 
and flap at two angles of attack, 0.47° and 5.65°, are shown in fig- 
ure 12(b). The predicted pressure distributions from both the deflected 
and undeflected wake methods presiime the gap between the wing and flap 
to be 0.02 c. The undeflected wake model predicts a slightly higher 
pressure difference on both wing and flap, confirming the previous 
results for lift coefficients, although the general shapes of the curves 
are very similar. 

A second set of results used for comparison purposes consists of a 
chordwise pressure distribution on a swept wing with a partial span flap 
(ref. 12) . This configuration, shown in the sketch in figure 13, has a 
35° leading-edge sweep angle, a flap span approximately one half the 
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semispan, and a flap chord equal to 20 percent of the wing chord. Calcu- 
lations were made at zero degrees angle of attack and 10° flap deflection 
angle. The lattice arrangement has 12 spanwise vortices by 8 chordwise 
on the wing and 4 chordwise on the flap. There is no gap between the wing 
and the flap. The data were taken at a Mach number of 0.6; but since the 
wing-flap program is incompressible, no correction was made for compres- 
sibility effects. The predicted pressure distribution from both vortex- 
lattice methods are in good agreement with each other as are the total 
lift and pitching-moment coefficients on the wing-flap configuration. 

The predicted pressure distribution has the same shape as the measurements, 
including the loading peak near the wing and flap junction. 

Finally, the effect of allowing the trailing legs of the horseshoe 
vortices to remain attached to the flaps is examined on a typical EBF 
transport model having a triple-slotted, full-span flap system (ref. 13) . 
The lattice arrangement used on this configuration consists of 15 span- 
wise vortices on the semispan and 4 chordwise vortices on the wing, 1 
chordwise on the first flap, and 2 chordwise vortices on each of the two 
aft flaps. The lattice elements are not evenly distributed in the span- 
wise direction, as this same lattice arrangement is used later for 
power-on calculations (fig. 26) . The model used in reference 13 has a 
large fuselage, four nacelles beneath the wings, and a leading-edge slat 
deflected 55°. None of these items have been accounted for in the 
predictions shown in figure 14. 

The measured and predicted lift and pitching-moment coefficients are 
shown in figure 14. As illustrated in figure 12 (a) , the effects of 
deflecting the trailing vorticity are to: (1) decrease the predicted 

lift coefficients over the entire range of angles of attack, (2) decrease 
the predicted lift curve slope a small amount, and (3) increase the nose- 
up pitching-moment coefficients. The net result of allowing the trailing 
vorticity to follow the flaps is to improve agreement with experimental 
results for both lift and pitching moments. 

The above comparisons between measured and predicted longitudinal 
aerodynamic characteristics of wings with flaps indicate that better 
accuracy is to be obtained by allowing the trailing vorticity to deflect 
with the flap system. In particular, the predicted pitching-moment 
coefficients and moment curve slope tend to be in better agreement with 
experiment when the deflected wake model is used, unless otherwise 
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specified, all further results presented in this report will correspond 
to the deflected wake prediction method. 

Jet Wake Model 

Single ring characteristics .- Since the complete jet is made up of 
a series of vortex rings, it is appropriate that the characteristics of 
a single, elliptic vortex ring be examined. The induced velocity field 
due to a circular vortex ring is well documented in reference 14; 
however, similar results for an elliptic vortex ring are not available 
so far as the authors can determine. In figure 15, the induced axial 
velocity in the plane of the ring is shown for a series of elliptic 
rings with varying axis ratios. The axial velocity distribution along 
both the minor (fig. 15(a)) and major (fig. 15(b)) axes are very similar 
and illustrate the singular behavior associated with induced velocities 
close to a vortex filament. The velocity distributions are similar for 
a wide range of axis ratios and there is effectively little change beyond 
an axis ratio of 10. The axial velocity outside the ring on the major 
axis tends to become very small as the axis ratio increases, because the 
ring configuration approaches two parallel vortex filaments with little 
periphery to induce a velocity outside the end of the ring. 

The induced axial and vertical velocities outside the plane of the 
vortex rings are easily calculated using the analysis presented in the 
appendices. Some induced velocities at various points on the vertical 
plane of symmetry are shown in Tables I and ll for elliptic vortex rings 
of different axis ratios. The limits of a circular ring and two infinite 
vortex filaments are also shown in the tables. There appear to be some 
numerical problems in calculating induced velocities in the plane of 
high-aspect-ratio elliptic vortex rings (a/b > 5) near the edge of the 
ring. There is a singularity on the ring itself, and calculated veloci- 
ties near the ring show some inconsistency. This is not significant, 
since this region of a vortex ring is never used in computing jet induced 
velocities. As described in reference 3, field points which fall between 
two adjacent rings but very close to one are treated as if they lie off 
the plane of the vortex ring exactly half way between two adjacent rings. 
This correction is carried out to achieve smooth induced velocity profiles 
which are in agreement with the profiles which would be obtained from a 
uniform distribution of vorticity on a semi-infinite cylinder. 
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Jet characteristics .- The jet model is made up of a series of coaxial 
vortex rings of elliptic shape, and the induced flow field at any point 
is determined by adding the contribution of each ring in the jet model. 
Even though the induced velocity from a single ring is highly nonuniform 
as shown in figure 15 and Tables I and II, the total induced velocity 
field associated with the complete jet is remarkably smooth and predict- 
able. The jet induced axial velocity profiles at various points inside 
a circular jet are discussed in reference 3 and illustrated in figure 5. 

It is seen that the axial velocity profile inside a circular jet is 
uniform inside the boundary of the jet and nearly zero outside the 
boundary. Several elliptic cross-section jets with various axis ratios 
are described as follows. In each of the jets to be described, the jet 
length and spacing between rings is the same, and the jet is allowed to 
expand to four times its initial area. 

A sketch of a typical jet model used in the following discussion is 
shown at the top of figure 16(a). The cross-section shape of each jet 
is similar over its entire length as the axis ratio is kept constant. 

In figure 16(a), predicted axial velocity profiles in the vertical plane 
of symmetry for three jets are shown. The jet cross sections vary from 
circular (a:b = 1:1) to moderately elliptical (a;b = 3:1), and the axis 
lengths are chosen such that the cross-sectional areas of the three jets 
are the same at all axial stations along the centerline. The profiles 
shown in figure 16(a) are computed at the end of a straight section of 
jet (Xj = 5 in the sketch) before any expansion has occurred. As was 
the case in reference 3, the profiles are nearly uniform across the jet, 
and the only difference between the circular jet and the two elliptic 
jets occurs near the edges of the jets. 

Predicted axial velocity profiles across a series of elliptic jets 
at a point where the local perimeter (P) is 1.25 times greater than the 
initial perimeter (Pq) are shown in figure 16(b). These profiles are 
also shown in the vertical plane of symmetry. The magnitude of the 
maximiim velocity in the jet is nearly the same for each of the various 
shape jets, but since some expansion has taken place, the magnitude 
is less than that shown in figure 16(a). Two very high-aspect-ratio 
elliptic jets, 5:1 and 10:1, are considered in figure 16(b). Their 
velocity profiles are similar to the other jets with the exception that 
the nonuniformity near the edges of the jets appears to increase as the 
cross sections become more elliptical. For all practical purposes, the 
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axial velocity profiles may be considered uniform inside the boundaries 
of the elliptic jets considered thus far. 

Predicted axial velocity profiles in the horizontal plane of symmetry 
of two elliptic jets are shown in figure 16(c). These profiles are 
similar to those just described except that they are more uniform across 
the width of the jet than across the height. The induced axial velocity 
tends to increase slightly from the centerline to the jet boundary and 
there is less difference between profiles as the jet aspect ratio 
increases . 

Finally, a series of profiles at various lateral stations across a 
3:1 elliptic jet are shown in figure 16(d). As would be expected, these 
profiles are all similar. 

The axial velocity profiles presented in figure 16 are typical of 
the results obtained from the vortex ring jet model, varying the expan- 
sion rates does not change the shape of the profiles, nor does changing 
the cross-sectional shape of the jet. For example, a jet may begin 
circular and gradually change to elliptical with distance downstream and 
the local profiles will still appear as those just discussed. 

As the jet expands with distance along the centerline, the induced 
axial velocity on the centerline decreases. The axial velocity is 
inversely proportional to the local perimeter of the jet, as illustrated 
in figure 17, where the axial velocity, normalized by the initial vortex 
cylinder strength, is shown as a function of the perimeter, normalized 
by the initial perimeter of the jet. A nxrnber of jets, ranging from 
circular to a 10:1 elliptical jet, are shown on this correlation curve, 
various expansion rates are represented by the points in this figure, 
from nonexpanding to very large expansion rates. There is some small 
deviation from the correlation curve as the jet expands to a perimeter 
approximately double its original perimeter. 

Using the information in figure 17, it can be shown that the mass 
inside the boundaries of an expanding circular vortex ring jet model is 
not constant but is increasing with distance along the jet as long as 
the perimeter of the jet is increasing. This is illustrated in figure 18 
where the predicted ratio of the local mass flow to the initial mass flow 
in the jet is shown versus the perimeter ratio. This curve was obtained 
from the same series of jets illustrated in figure 17. Figure 18 indi- 
cates that at any point along an expanding jet, the mass ratio is deter- 
mined by the local perimeter of the vortex rings; or conversely, if a 
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specific mass ratio at a point in a jet is desired, it can be obtained 
by choosing the correct expansion of the vortex rings. 

Measured axial velocity profiles in the wake of a JT15D-1 jet engine 
mounted beneath a wing are available in reference 13 for a range of thrust 
levels. Two profiles corresponding to thrusts of 5338 newtons (1200 lbs) 
and 3114 newtons (700 lbs) are shown in figures 19(a) and (b) , respectively. 
The profiles were measured on both the upper (wing) side and the lower 
(free) side of the engine centerline (assumed straight) at a point approxi- 
mately two wake diameters downstream of the engine exit. The data indicate 
that the wake is nearly symmetric about the centerline; therefore, one half 
of the profiles is illustrated with a faired curve. Since the velocity 
profile was measured in a single vertical plane through the centerline, 
we must assume the jet to be axisymmetric . 

A circular jet model was designed to expand at a rate that would 
produce the same mass flow at the measured profile station, and the 
resulting predicted velocity profile is shown dashed in figures 19(a) 
and (b) . A comparison of the momentum in the measured and predicted jets 
indicates that the real jet has approximately 5 percent more momentum 
than the predicted jet. Therefore, it is possible to match satisfactorily 
both the mass and momentum of the real jet with a vortex ring model. Any 
interference calculations taking place inside the jet boundaries will be 
reasonably accurate when averaged over the area of the complete wake, but 
there may be certain inaccuracies locally inside the wake due to diff- 
erences in the shape of the velocity profile. This result will be 
illustrated in more detail when the predicted span loadings on an EBF 
configuration are compared with data. 

Centerline calculation .- An important feature of the EBF prediction 
method is the calculation of the deflection of a jet due to the influence 
of the flow field induced by the loading on the wing and flap surfaces. 
Ideally, the centerline calculation procedure described previously 
should be checked out by comparing the predicted path of a jet in a 
nonuniform velocity field with the actual path, but data of this type 
are not available. 

A large amount of jet path data are available in reference 15 for 
different strength jets directed into a uniform crossflow at various 
angles. These data are summarized in the following empirically 
determined equation: 
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which is equation (4) of reference 15. 

The centerline calculation method was applied to this jet flow case 
to evaluate the iteration procedure. The vortex ring model of the jet 
was based on a typical axisyrametric jet expansion rate and the jet 
centerline was assumed to be straight initially. Predicted centerlines 
for two jets directed at a 30° angle to the free stream are shown in 
figures 20(a) and (b) for jet velocity ratios 8.85 and 2.81, respectively. 
In both cases, the predicted centerlines agree very well with the center- 
lines given by equation (41) . 

The agreement starts to deteriorate beyond ten jet diameters down- 
stream of the origin. This is approximately the end of the initial 
region of a jet where the potential core disappears. The jet expansion 
rate changes at this point, which changes the jet velocity on the center- 
line. The lack of agreement downstream of this point is probably related 
to some inaccuracy in specifying the expansion rate of the vortex ring 
jet model. 


Convergence Characteristics 

The previous sections have been concerned with the individual nature 
and accuracy of the wing-flap and jet wake methods. Since the EBF pre- 
diction method is a combination of these two methods with iteration to 
obtain the mutual interference, it is of interest to examine the 
convergence characteristics for an EBF configuration. For this purpose, 
calculations were made on the four-engine configuration of references 13 
and 16, which has a swept wing with a triple-slotted full-span flap 
system. The calculations were made with the computer program of refer- 
ence 4 using a given lattice arrangement and a circular cross section 
expanding jet (described in the following section) . The initial assump- 
tion for the jet centerline in all cases is a straight line coincident 
with the engine centerline. Convergence characteristics were examined 
for total normal force, individual flap normal force, span-load 
distribution, and centerline position. 

Convergence of the procedure for wing-flap normal-force coefficient 
is illustrated in figure 21. The normal-force coefficient from the first 
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iteration is approximately 50 percent larger than the results from the 
second iteration. After four full iterations, both curves, representing 
thrust coefficients of 2.3 and 4.0, have converged to within 7 percent. 
Additional iterations would produce more nearly converged results, but 
7 percent was judged low enough to demonstrate convergence. The conver- 
gence pattern for normal-force coefficient is the same for this configu- 
ration at other angles of attack and for other similar configurations. 

The convergence of the normal-force coefficient on each lifting 
surface of the wing-flap configuration is shown in figure 22. In this 
case, the normal force is defined to be normal to each individual compo- 
nent. Each component tends to converge according to its own pattern, but 
all components reach convergence at about the same rate. Flaps 2 and 3 
have nearly identical convergence patterns because they have similar 
orientation to the jet wakes. The wing and flap 1, both with zero inci- 
dence, show a small effect of wake centerline deflection on their normal 
force because they are aligned parallel to the initial centerlines and 
small vertical motions of the centerlines have little effect on the 
induced velocity fields on these surfaces. 

The convergence pattern of the spanwise distribution of normal force 
on each lifting surface is shown in figure 23. The wing normal-force 
distributions through four iterations are presented in figure 23(a). The 
peaks in the loading near t] = 0.25 and 0.42 are caused by the engine 
wakes moving upward due to angle of attack and interfering directly on 
the aft portions of the wing. The normal-force distribution is nearly 
converged after two iterations as there are only small corrections pro- 
duced by iterations three and four. The maximum normal-force coefficients 
at the peaks show the largest change between successive iterations, but 
these large loadings act on such a small portion of the wing that the 
total wing normal force is only slightly affected by the peak loadings. 

Corresponding normal-force distributions on flaps 1, 2, and 3 are 
shown in figures 23(b), (c) , and (d) , respectively, for the four iterations 

Flap 1 has high peak loadings similar to those on the wing which is not 
unexpected since flap 1 has zero deflection and can be looked at as simply 
an extension of the wing. The high peaks corresponding to the first 
iteration are caused by the edge of the jets hitting one row of control 
points. As the jet is bent down by the iteration procedure, the loading 
is reduced by the third iteration and has obviously converged after the 
fourth iteration with the peak loadings reduced considerably from their 
initial values. 
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Flap 2, with a deflection angle of 20°, is deflected into the center 
of the initially straight jet. Thus, the jet wakes are directly 
influencing a larger area of this flap than on the previous flap. This 
is illustrated by the width of the peak loading regions on flap 2. Each 
peak covers approximately 10 percent of the semispan of the flap. With 
each iteration, the jet is bent downward to approach the flap angle and 
the peak loadings are reduced, but the area of influence of the jets 
remains larger than on flap 1. As shown in figure 22, the total loading 
on flap 2 is much larger than that on flap 1. 

The normal-force distribution on flap 3 after one iteration is quite 
different in character than the distribution observed on the other flaps. 
This loading has a single peak which is higher than that on any other flap. 
The single high peak is caused by the overlapping of the straight, expanding 
jets. The flap is deflected 40° into the center of the overlap region of 
the wakes where the total induced effects are nearly double those on 
preceding flaps, thus resulting in inordinately high local loadings. With 
iteration, the wakes are deflected downward such that the overlap region 
occurs beneath the flaps. This can be seen in the third and fourth itera- 
tions where the twin peak character of the flap loading has returned. The 
converged normal-force distribution on flap 3 has peaks which are broader 
but lower in maximum magnitude than those of flap 2 because the jets have 
expanded to a larger radius with the resultant lower jet velocities. 

Thus, flap 3 experiences lower jet induced velocities over a larger portion 
of the flap than the flaps upstream. 

The convergence of the jet centerlines is shown in figures 24(a) 
and (b) for the inboard and outboard jets respectively. The initial jet 
boundary is also shown in these figures as it is superimposed on the 
wing-flap model. Note that the boundary is straight, indicating a linear 
expansion from the engine exit. The initial centerline used for the 
first iteration is straight as shown. The combination of the jet induced 
velocities, the wing-flap induced velocities, and the free stream produce 
the centerline, labeled as number 2, which is used for the second 
iteration. As was inferred from the loading distributions, convergence 
of the centerlines has been achieved by the fourth iteration. The final 
centerline calculation at the end of iteration 4 is denoted as centerline 
niamber 5 which, in each figure, is nearly the same as number 4. Thus, 
the fifth centerline is the converged centerline for each jet. The jet 
boundary can be placed on this centerline and the total flow field 
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calculated. Such a flow field is compared with a measured flow field in 
the next section concerning comparisons with data. 

Lift and pitching-moment coefficients at two angles of attack, 
a = 6.5° and 18.5°, are shown in figure 25 for each iteration to illus- 
trate convergence of lift and pitching -moment curve. The values of lift 
coefficient used on these curves are for the wing and flap lifting 
surfaces only and do not include any lift components due to engine thrust, 
fuselage lift, or nacelle lift. Iteration 1 at each angle of attack was 
started with the same straight centerline jet model. The lift curve from 
the first iteration is nearly flat, indicating a very small effect of 
angle of attack. It is apparent that the jet induced effects overshadow 
all others when the centerlines are straight. At the end of the second 
iteration, the lift at a = 6.5° is less than that at a = 18.5°. Iteration 

2 also produces the greatest lift curve slope. The jet induced effects at 
the lower angle were reduced because the jet centerlines moved downward 
further from the lifting surface than did the same centerlines at the 
higher angle. This occurs because the higher angle of attack flow tends 
to keep the jet closer to the wing and flaps. Iterations 3 and 4 produce 
the same lift and moment curve slopes and the lift values at both angles 
of attack agree within approximately 7 percent. The converged jet center- 
lines calculated at a = 6.5° were used for the initial centerlines for the 
a = 0° calculation. In two iterations the solution had converged within 

3 percent. The converged result at a = 0° is also shown in figure 25. 

The convergence results described above are typical of those 
observed on other EBF configurations. The method has never failed to 
converge, although convergence is slower for high flap angles. Generally, 
calculations have been initiated with a straight centerline because of 
simplicity, but obviously the number of iterations would be reduced if 
the initial centerline more closely approximated the final one. On the 
basis of cases run, the convergence is more rapid if the centerline 
approaches its converged position from above rather than below, because 
the correcting velocities causing the centerline position to change are 
larger if the centerline is too close to the wing-flap. 

EBF Data 

The overall EBF prediction method was evaluated by comparing 
predicted results with data on several typical EBF configurations. 

Ideally, one would like to have wing-flap-engine data on overall forces 
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and moments, loads on individual flap segments, pressure distributions, 
and jet wake velocity surveys. A large amount of data is available, but 
most of it is on complete configurations and consists only of overall 
force and moment coefficients. Thus, there are generally uncertainties 
due to force contributions from the fuselage, tail, nacelles, and 
leading-edge slats and to lack of direct knowledge of engine thrust. The 
results presented in this report attempt to account for these additional 
factors in the . following manner. 

Corrections . - The predicted aerodynamic characteristics obtained 
from the EBF method described in this report are for a wing-flap- jet 
configuration only. No fuselage, nacelle, or leading-edge slat effects 
are included. Results in reference 16 indicated that these effects are 
not negligible and should be included with the predictions to give 
meaningful comparisons with the data. Thus, the following corrections 
are made to the predictions. 

The total lift coefficient is 

+ C, sin a + C_ + (42) 

L |J. J-lg Ijg 


where the first term is from the wing-flap vortex-lattice method with 
jet interference, and the second term is the component of total thrust 
in the lift direction. The last two terms in equation (42) are the 
combined lift of the leading-edge slat, the fuselage, and the nacelles. 
Typical values of these terms for the configuration under consideration 
are presented in Table III as a function of angle of attack and thrust 
coefficient. These were obtained in reference 16 from the difference 
between the total measured lift coefficients and the lift coefficients 
obtained from integration of the measured pressure distributions. 

The total pitching-moment coefficient is 
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where is the result from the vortex-lattice prediction method. 

WF 

The value is the pitching moment due to the force on the leading- 

O 

edge slat and is estimated by integrating the section normal-force 
distributions in reference 16. These results are also shown in 
Table III. The value the pitching moment due to the fuselage 

and nacelles, is estimated using slender-body theory to be approximately 
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Although the body contribution to pitching moment is small in magnitude 
compared to the wing-flap contribution, it can have an appreciable effect 
on the total moment curve. The moment center is such that there is no 
moment produced by the thrust of the engines. 

The total (nonviscous) drag coefficient produced by the prediction 
method is 
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where is the induced drag coefficient calculated by the vortex- 

WF „ 

lattice method on the wing and flaps. The value is the ram drag 

coefficient for the engines. This contribution is given in figure 4 of 

reference 13 for the particular configuration being considered. These 

results are reproduced in Table IV for convenience. The last term is 

the component of total thrust opposing the drag of the model and is 

included in the predictions to put them in the same form as the data. No 

estimate of viscous drag is included in the prediction method. 

Data available .- The best data available for comparison purposes are 
that of references 13 and 16. Gross forces and moments and pressure 
distributions were obtained on a four-engine EBF configuration; thus, the 
individual loading distributions on each lifting surface are available. 
Because these data are unique, they will be used almost exclusively to 
evaluate the prediction method and the associated computer program. 

The model is a large-scale, four-engine configuration with an aspect 
ratio 7.28 wing having a 25° swept quarter chord and a 0.4 taper ratio. 
High lift devices include full-span, triple-slotted, trailing-edge flaps 
and a leading-edge slat deflected 50°. Two flap deflections are consid- 
ered for the full-span flap configuration: a take-off setting of 0°, 

20°, and 40° for flaps 1, 2, and 3, respectively, and a landing setting 
of 15°, 35°, and 55°. The JT15D-1 engines are housed in nacelles mounted 
beneath the wing on pylons at the t] = 0.25 and 0.42 semispan stations. 
Tests were run at gross thrust coefficients between 0 and 4.0. More 
detailed configuration data are available in reference 13. 

Chordwise static pressure distributions at ten spanwise stations on 
each lifting surface are described in reference 16. Section normal-force 
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coefficients calculated from integrated pressure data on each lifting 
surface are presented in the same reference. These integrated pressure 
data are used for the following comparisons . 

The vortex-lattice arrangement for the wing and flaps is shown in 
figure 26 where the flaps are shown undeflected so that the relative size 
of all the lattice elements can be compared. The lattice layout is based 
on the guidelines presented in Appendix A of reference 3, where possible. 
One change that is required by the deflected trail ing-leg vortex-lattice 
method is that an equal number of spanwise elements be placed on each 
lifting surface and that all trailing legs be aligned with each other. 

The lattice shown in figure 26 is tailored for the two jet wakes at 
T] = 0.25 and 0.42 and should be considered a "minimum" lattice for a 
four-engine wing-flap configuration of this complexity. The effect of 
changing the lattice layout is investigated later in this report. 

A second set of data considered for comparison purposes is that for 
the four-engine EBF model of reference 17. This model has an aspect 
ratio 7, unswept, cambered wing with double-slotted, full-span trailing- 
edge flaps an<3 a leading-edge slat deflected 55°. The model engines are 
cold-gas-turbine-driven fans which simulate the wake of high-bypass-ratio 
turbofan engines. The engines were mounted at 3° incidence to the wing 
root chord in order to direct their exhaust upward toward the lifting 
surfaces (position 4 in ref. 17) . The leading-edge slat has a chord 
equal to 19 percent of the reference chord and is treated as a highly 
cambered portion of the wing. The two flaps have chords of 15 and 30 
percent of the reference chord. This is the same configuration considered 
in the comparisons of reference 3 . 

Circular Jet Results 

The jet wake model used in the initial calculations has a circular 
cross section over its entire length. The vortex cylinder strength, 

Y/V, is obtained from equations (28) and (29) , and the expansion rate 
is determined from figure 11 which fixes the jet cross-sectional area 
aft of the last flap. A linear variation of wake radius is used to 
specify the jet growth between the engine exit position and the axial 
station defined by figure 11. The jet radius is extrapolated at the 
same expansion rate over the remainder of the length of the jet. 

Take-off configuration .- The first configuration on which comparisons 
with data are made is the four-engine model with take-off flap setting 
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(6£ = 0°/20°/40°) described in references 13 and 16. The jet turning 
efficiency was assimed to be 85 percent which limited the jet deflection 
angle to -34° for all calculations. The predictions to follow have all 
converged to within approximately 8-percent tolerance. The convergence 
is not the same at all angles of attack, which introduces a small 
uncertainty in the slope of the predicted curves. 

In figure 27 the predicted section normal-force coefficients on 
flaps 1, 2, and 3 are compared with experimental results at C = 4 and 
a = 18.5° obtained from reference 16. Wing data are not available for 
this configuration. The predicted values are from the fourth iteration 
described in figure 23. Figure 27(a) illustrates that the predicted peak 
loadings due to jet wake interference on flap 1 are greater than those 
measured and cover a smaller portion of the wing. The dip in the pre- 
dicted loading between the jets is less than the measured values as is 
the predicted loading outboard of the second wake (t) > 0.45) . This 
indicates that the model jet has not expanded sufficiently at this 
station and perhaps should be expanded to a larger cross-sectional area. 

As noted on the figure, the predicted total normal-force coefficient on 
the single flap is approximately 30 percent higher than the value 
obtained by integrating the measured distribution. 

The comparison in figure 27 (b) for flap 2 shows the same pattern in 
that the predicted normal force in the vicinity of the engine wakes is 
much larger than the measured quantities with the result that the total 
force on the flap is overpredicted by nearly 40 percent. 

The comparison for flap 3 in figure 27 (c) shows good agreement 
betv/een the predicted and measured loading distributions. The peak 
loadings, the width of the loading, and the total force on the flap are 
all in good agreement. Since this flap is nearest to the point at which 
the jet wake area is specified, based on measurements near a similar EBF 
configuration, the jet model is probably in better agreement with the 
actual jet on this flap than on the previous two flaps. 

The predicted and measured longitudinal, aerodynamic characteristics 
on the four-engine EBF model with take-off flap configuration are compared 
in figures 28(a) and (b) . The predicted lift, pitching -moment, and drag 
results are obtained from equations (42) , (43) , and (45) , respectively. 

The predicted power-off lift curve has the correct slope but is about 
10 percent higher than the data. The corresponding pitching-moment curve 
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also has the correct slope, but the predicted results have a larger nose- 
down moment than the data. The predicted drag coefficient in figure 28(b) 
is slightly lower than the data as would be expected since viscous drag 
is not considered in the prediction method. 

The lift coefficient curve in figure 28(a) for power-on, C = 2.32, 
is lower than the data, but the predicted slope is higher. It appears 
that the method is converging on a lift coefficient that is too low at 
low angles of attack. This result may be caused by the converged jet 
being too far beneath the lifting surfaces, a difficilty caused by 
estimating too high a value for jet turning efficiency. The predicted 
pitching-moment curve generally has too little nose-down moment in the 
region of low lift, indicating that the lost lift is aft of the moment 
center and most likely the lost lift should be on the flap surfaces. The 
drag curve shows very good agreement with the experimental results . 

The predicted results in figure 28(a) corresponding to the highest 
thrust coefficient, = 3.95, indicate the correct level of lift but a 
lift curve slope slightly higher than that measured. As the case for 
the lower thrust coefficient, the region which has the best agreement 
between measured and predicted lift coefficients also has best agreement 
for pitching-moment results. The drag coefficients are in fair agreement 
with the data. 

Landing configuration .- The next set of experimental results used 
for comparison with predictions is on the same EBF model in a landing 
configuration with flaps 1, 2, and 3, deflected 15°, 35°, and 55°, 
respectively. In addition to overall force and moment coefficients, 
data are lavailable on span— load distribution and chordwise pressure 
distributions from reference 16. The same vortex-lattice arrangement 
(shown in figure 26) was used for these calculations. As shown in ref- 
erence 13, the static turning effectiveness of the jet acting on this 
configuration is approximately 0.70; therefore, the jet turning angle 
is limited to a maximum deflection of -38.5°. 

The results to follow are all converged within approximately 11 
percent. Most of these results are obtained from the fourth iteration 
of the solution, the same niimber of iterations used for the take-off 
configuration results. As discussed in those results, not all angles 
of attack are converged to the same tolerance; therefore, some uncertainty 
as to the slope of the predicted curves exists . 
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In figure 29, measured and predicted spanwise distribution of 
section normal-force coefficients are compared for the landing-flap 
configuration at 18.5° angle of attack. Two power settings, power off 
and C =2.3, are shown in this figure. Results on the wing, shown in 
figure 29(a), indicate that the predicted loading for power off is higher 
than was measured to the extent that the total predicted normal force is 
almost doiible the measured value. When the engines are operating, the 
predicted wing loading is in good agreement with experiment in the near 
vicinity of the engines; however, outboard of the engines (t] > 0.5) the 
predicted results are too low. Integrating the two curves gives a 
predicted normal-force coefficient, on one wing panel, of 1.37 to be 
compared with the experimental result of 1.57. 

Similar results on flap 1 are shown in figure 29(b). The predicted 
pov/er-off loading distribution is again above that measured with the 
result that the total predicted normal force is much larger than the 
measured result. With power on, the predicted interference loading due 
to the inboard engine wake is higher than the measured value and covers 
a more narrow band of the wing than is illustrated by the measurements. 
The opposite situation is evident near the outboard engine; and outboard 
of the engines, the predicted loading is less than the measured results. 
The difference in levels of the predicted peak loadings near the two 
engine centerlines is due to small differences in the inboard and out- 
board jet centerline heights with respect to flap 1. This. causes the 
jet induced velocities at the control points in line with the two center- 
lines to be different, thus producing quite different loadings. 

Measured and predicted loading distributions on flaps 2 and 3 
are shown in figure 29(c) and (d) , respectively. The power-off results 
are similar to those on the previous flap. The predicted power-on 
loading curves on both flaps exhibit similar shapes, and on each flap, 
the loading around the engine centerlines is similar in shape and 
magnitude. The measured loadings on each flap show a larger interference 
effect from the outboard engine than from the inboard engine. The 
predicted total loading on each flap is in reasonable agreement with 
the measured result. As is the case on all flap surfaces, the predicted 
loading drops off more rapidly outboard of the mid-semispan than the 
measurements indicate. 
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Measured and predicted section normal-force coefficients on the 
same configuration at 18.5° angle of attack are compared in figure 30 
for a high thrust setting, C = 4. The only difference in the theoretical 
model used for these results and the model used for the lower thrust 
setting in the previous figure is the jet model. The higher thrust jet 
has a higher jet exit velocity; consequently, the jet expands at a 
faster rate than the jet corresponding to C = 2.3. 

The effect of the change in jet model is seen on the predicted 
results on the wing shown in figure 30(a). The predicted span loading 
on the wing has two loading peaks which are caused by the direct inter- 
ference of the jet on the wing lifting-surface model. Similar peaks 
are not shown by the data. The differences between the measured and 
predicted loadings on the wing in the vicinity of the jet wakes will 
be discussed later in this section when chordwise pressure distributions 
are examined. The predicted loading on the outboard portion of the wing 
is again lower than measured results; however, the additional loading 
due to the peaks causes the predicted net normal force on the wing to 
be in good agreement with the value obtained from integration of the 
experimental curve . 

The normal-force distributions on flap 1 are shown in figure 30 (b) 
where the predicted peak loadings are again much larger than those 
measured. The peak attributed to the outboard jet is at » 34 which 

is off the scale chosen for this figure. 

Results on flap 2, figure 30(c), show good agreement for the loading 
induced by both the inboard and outboard jets, except that again the 
predicted loads drop more rapidly outboard of the outer jet than is 
indicated by the data. The total normal force on this flap is in good 
agreement with the value obtained from integration of the experimental 
curve . 

The measured and predicted normal-force distributions on flap 3 are 
compared in figure 30 (d) . The measured results appear as on previous 
flaps, but the predicted results have a different appearance brought 
about by the lack of two well-defined areas of peak loading. This is 
caused by the meeting of the jet boundaries in the vicinity of the third 
flap. The expansion rates are such that the two jets come together near 
this flap causing the jet induced loading to be spread out over the 
entire interference region of the flap. The small dip occurring at 
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= 0.37 is due to one of the two control points at this spanwise station 
falling just outside the boundary of the outboard jet, thus reducing the 
induced velocity field at this point. 

Similar comparisons of the measured and predicted section normal- 
force coefficients on the EBF configuration of references 13 and 16 with 
the flaps deflected 15°, 35°, and 55° and a high thrust condition are 
presented in figure 31 for a lower angle of attack (a = 6.5°) . The jet 
model is identical in strength and expansion rate to the model used for 
the previous figure, but the jet centerlines converge to different posi- 
tions due to the change in angle of attack. The results have the same 
overall appearance as those in figure 30. 

The chordwise distribution of pressure coefficients on the wing and 
flap surfaces at a semispan station aligned with the outboard engine are 
shown in figure 32. The data are shown as dashed curves with the long 
dashed portion of each curve representing the length over which measure- 
ments were made and the short dashed portion representing extrapolation 
to the trailing edge. Since the vortex-lattice method implies constant 
pressure on each lattice element, the theoretical pressure coefficients 
are shown as such. Too few chordwise lattice elements are used to permit 
a smooth cruve through predicted points. Notice that the jet induced 
loading on the aft quarter of the wing chord contributes all the addi- 
tional loading Which makes up the peak shown in figure 31(a) . The 
absence of pressure taps on this part of the wing does not permit evalua- 
tion of the predicted pressure rise. The airfoil section in this region 
of the wing is distorted to provide for retraction of the flaps; therefore, 
it is unlikely that an idealized lifting-surface theory with no thickness 
effects included can adequately predict the detailed loading distribution 
in this region. 

The measured and predicted longitudinal aerodynamic characteristics 
of the four-engine EBF model with landing flap configuration are compared 
in figure 33(a) and (b) . The predicted results presented on this figure 
are corrected for the body, slat, and nacelle effects using the same 
information used to correct the take-off configuration results. These 
results are similar to those shown for the lower flap deflection angles 
in figure 28. The predicted power-off lift curve has the correct slope 
but is above the data. The power-on results are within 10 percent of 
the data over the angle of attack range shown (0 < a < 18.5°) . 
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The pitching-moment results shown in figure 33 agree reasonably 
well. The inconsistent shape of the predicted curve at C =2.3 may 
be due to differences in the degree of convergence of the calculations 
at different angles of attack, since the loadirig, particularly the 
pitching moment, is sensitive to small changes in wake position. 

The drag results are shown in figure 33(b). The results are 
generally similar to those at the lower flap angle. One would expect 
best agreement at the higher thrust settings and negative drag values 
(small angles of attack) , since the power effects are largest relative 
to viscous effects at these conditions. Such is the case, both with 
figures 33(b) and 28(b). 

Flow fields .- The predicted flow field aft of the trailing edge of 
the last flap at a spanwise station corresponding to the centerline of 
the inboard jet is shown in figure 34. In figure 34(a), the final con- 
verged flow field at a = 6.5° and C = 3.95 for the take-off flap 

M' 

configuration (6^ = 0°/20°/40°) is shown. The uniformity of the jet flow 
characteristics of the vortex ring model is well illustrated. In the 
inset, a measured flow field aft of a similar EBF configuration under 
similar flow conditions is reproduced from reference 7. The measured 
flow field, also aligned with the centerline of the inboard jet, is very 
similar in magnitude and direction to the predicted flow field. The 
same type of comparison and results are shown in figure 34(b) for a 
lower thrust coefficient and higher angle of attack. 

Lattice arrangement .- The effect of increasing the number of lattice 
elements representing the wing-flap configuration of references 13 and 16 
has been investigated. Additional chordwise rows of vortices were added 
to each flap (fig. 26) so that flap 1 had two rows, and flaps 2 and 3 
each had three rows of vortices. Results from a power-on calculation 
indicated less than 2-percent difference in the predicted aerodynamic 
coefficients. As anticipated from the conclusions of reference 3, 
increasing the number of lattice elements reduced the magnitudes of the 
predicted forces. Based on this investigation, it is felt that the 
lattice shown in figure 26 is probably a reasonable lattice for the 
particular configuration for force calculations. If more detailed 
pressure distributions are desired, one would want to increase the 
number of chordwise area elements. Guidelines for choosing a lattice 
layout for EBF configurations are presented in the program manual 
(ref. 4) . 
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Unswept wing configuration .- The prediction method is next applied 
to the rectangular wing, four-engine, EBF model of reference 17. 
comparisons of measured and predicted longitudinal aerodynamic character- 
istics on this model with trailing-edge flaps deflected 17.5° and 35° 
are shown in figure 35(a). The power-off predicted lift ciurve is approxi- 
mately 10 percent above the data, but it has the correct slope. The 
moment curve also has the correct slope, but too much nose-down moment 
is predicted. The lift and moment curves are corrected for the force 
and moment on the fuselage using slender-body theory. The predicted drag 
curve, v;hich includes no friction drag, is too low as would be expected. 

The power-on comparisons shown in the same figure are for a gross 
thrust coefficient of 5.5. For these predictions, the jet deflection 
angle limit was taken as -22° from the static turning curves in refer- 
ence 17. The comparisons between the measured and predicted results 
shov/ excellent agreement for both lift and pitching-moment curves, but 
the predicted drag curve is again lower than the data. The predicted 
lift and moment curves are corrected for the presence of the fuselage 
v;ith the same factors used in the power-off case. The lift and drag 
curves are also corrected for thrust effects, and an estimate of the 
engine ram drag coefficient is included in the predicted drag curve. 

The same configuration with flaps deflected 27.5° and 55° is shovm 

in figure 35(b).l The power-off lift is overpredicted by approximately 

30 percent, but ^he lift curve slope is in good agreement with the\ experi- 

u, 

mental value. The only explanation for the lift results is that there 
may be some separation on the model with the large flap deflection. The 
measured moment curve indicates a larger nose-up moment component than 
the prediction. This is indicative of a loss of lift on the flaps. 

When the engines are operating at a gross thrust coefficient of 
5.5, the predicted results are approximately 10 percent beneath the data 
at low angles of attack. The predicted lift curve slope is also higher 
than the measured slope. The pitching moments are in reasonably good 
agreement. The predicted results for the higher flap deflection were 
obtained with the same engine wake model used for the previous figure, 
but in this case, the jet deflection angle was limited to approximately 
-38° as recommended by the static turning curves in reference 17. This 
may be an overly optimistic turning efficiency, and a lower choice would 
have the effect of raising the predicted lift curve in figure 35(b). 
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Elliptical Jet Results 

The results presented thus far have all been obtained using a 
circular cross-sectional jet model. The effect of using an elliptical 
cross-section jet will now be examined on the four-engine EBF configu- 
ration with landing flap setting of references 13 and 16. The elliptic 
jet used for these calculations has the same momentiim and the cross- 
sectional area distribution along the jet centerline as the circular jet. 
The elliptical jet is assxmed to have an axis ratio of 2:1 at a point 
just aft of the trailing edge of the last flap. Since no data are 
available on the manner in which the jet axes change from circular to 
elliptical, a linear variation is chosen. The details of the specific 
elliptic cross-section jet used in the following calculations are 
presented in the program manual (ref. 4) . In an effort to conserve 
computation time, the initial guess for the elliptic jet centerline is 
taken as the circular jet centerline two iterations before convergence. 
The elliptic jet calculation was carried out for two iterations; and all 
comparisons which follow are taken from the last iteration of both the 
circular jet and elliptic jet results. The same jet turning efficiency 
used in the previous circular jet calculations was used for the elliptic 
jet calculations. 

The predicted section normal-force distributions corresponding to 
the circular jet and elliptic jet results are compared with experimental 
distributions in figure 36. The circular jet curves and the data are 
the same information shown in figure 30. The predicted loading on the 
wing under the influence of the elliptic jet is shown in figure 36(a). 

The interference loading due to the elliptic jet is less than the loading 
due to the circular jet. The major difference is the absence of the peak 
loadings which are caused by the circular jet moving up to contact the 
wing. The elliptic jet, in the vicinity of the wing trailing edge where 
the peak loadings were concentrated, has a smaller vertical height than 
the circular jet and it does not interfere directly on the wing. As 
shown on the figure, the total normal force on the wing under the 
influence of the elliptic jet is less than the measured result. 

Similar loading comparisons on flap 1 are shown in figure 36(b). 
Again the elliptic jet interference is less than the circular jet results, 
and there is a large difference in the predicted distributions in the 
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region of the outboard jet (t) == 0.42). The explanation for the differ- 
ences between the two predicted loadings is that the elliptic jet, being 
of different shape, nearly misses the control points near t] = 0.42 on 
this flap while the circular jet has a shape which results in direct 
interference on control points in this region. As was the case on the 
wing, the total normal force on flap 1 under the influence of the 
elliptic jet is less than the measured value. 

The circular and elliptic jet results have a similar shape on flap 2 
as shown in figure 36(c). The elliptic jet induced loading is less than 
the circular jet loading in the region near t] = 0.30. This is again 
caused by the elliptic jet missing certain control points. However, the 
total results from both jets are in agreement with each other and the 
data. 

An unusual predicted loading distribution is produced on flap 3 
by the elliptic cross-section jet. The high peak loading at t] ~ 0.34 
is caused by the overlapping of the inboard and outboard jets directly 
at one control point. This causes the loading on one area element to 
be double its usual level, thus resulting in an unrealistic loading 
distribution and a large net loading on this flap. 

The above results indicate some large differences in local loadings 
due to a change in jet cross-sectional shape. Some of these differences 
are caused by the choice of the vortex-lattice arrangement. The same 
lattice was used for both calculations, and it is possible that the 
lattice could be optimized for an elliptic jet shape. For example, on 
flaps 1 and 2, several control points are just outside the boundary of 
the jet model; consequently, the entire lattice element is unaffected 
by the direct interference of the jet. 

It is interesting to compare the total normal force on the complete 
wing-flap configuration for each jet shape. The normal -force coefficient 
on the circular jet configuration is 7.85, and that on the elliptic jet 
configuration is 7.69, a difference of only 2 percent. The loading is 
distributed differently over the wing and flap surfaces, but the total 
normal force imparted to the wing-flap configuration by the two jets is 
nearly the same. Since both the circular and elliptic jets had the 
same momentum and assumed turning efficiency, it can be concluded that 
similar total lifting-surface loading is required to turn the jets 
beneath the wing and flaps. The cross-sectional shape of the jet is 
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important if local loading and pressure distributions are required, but 
if gross aerodynamic forces are the goal of the calculation, the jet 
cross-sectional shape is not as important as the momentum and the 
entrainment of the jet model. 

CONCLUDING REMARKS 


An engineering prediction method has been developed to predict the 
loading distributions and longitudinal aerodynamic characteristics of 
externally blown flap configurations. Two potential flow models are 
combined for the method. The first is a vortex-lattice lifting-surface 
method capable of representing a wing of arbitrary sweep, taper, dihedral, 
camber, and twist with large multiple-slotted, partial or full-span flaps 
which can have large deflection angles. The second is a vortex ring model 
of a jet wake from a turbofan engine which represents the mass, momentum, 
and spreading of a turbulent jet with circular or elliptic cross section. 
It may be placed along a curved centerline. These models are combined in 
an iterative fashion to predict the loading on the wing and flaps in the 
presence of the jet wakes. The method yields chordwise pressure distri- 
butions and span-load distributions on each individual lifting surface 
and the total forces and moments on each surface. The method has the 
capability to compute flap hinge moments, but no results are presented 
in this report because of the lack of data for comparison purposes. Use 
of the method is restricted to the case of symmetry of the configuration 
about the vertical plane through the wing root chord. 

The complete prediction method was applied to EBF configurations on 
which detailed aerodynamic data are available. The comparisons of 
measured and predicted gross lift, drag, and pitching-moment coefficients 
on configurations with moderate flap angles (6^ < 40°) indicate generally 
good agreement for all thrust levels. This is felt to be due principally 
to the correct modeling of the entrainment and momentum characteristics 
of the engine wakes and to the proper treatment of the mutual interfer- 
ence between the jet wake and wing-flap. The interference model creates, 
on the wing-flap, both the momentum reaction due to jet deflection and 
the additional induced circulation characteristic of externally blown 
flap systems. 

As the flap angles increase beyond 40°, the predicted results agree 
less well with the data for both swept and unswept wings. As was 
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indicated by the intermediate iteration solutions, the predicted loads 
are sensitive to centerline location. The assumption that the wing-flap 
induced interference on the jet affects only its centerline and not its 
boundary becomes less accurate as the jet is more highly deformed, and 
it is probable that this is responsible for the poorer agreement at the 
higher flap angles. 

Comparisons of measured and predicted spanwise distributions of 
section normal force on each lifting surface indicate good quantitative 
agreement in some cases and poor agreement in other cases. Generally, 
the correct qualitative behavior is shown in which large peak loadings 
occur locally on the flaps due to direct impingement of the jet wakes, 
but the magnitude of the peaks is not consistently in agreement with the 
data. The differences are felt to be due to a combination of three 
factors: vortex-lattice arrangement, modeling of the velocity profile 

within the jet wake, and the boundary of the wake. 

The distribution of loading on the wing and flap segments is 
sensitive to lattice layout, in that a control point may be just inside 
or outside the wake impingement region and have quite a different load 
on its panel than that of an adjacent panel. As the number of panels 
is increased, sensitivity of the total loading to the loading on indivi- 
dual panels decreases; therefore, the complete solution becomes less 
sensitive to lattice/ jet location characteristics. 

The jet wake model has the proper mass, momenttim, and induced velocity 
field outside the jet, but a velocity profile inside the jet which is 
probably quite different than the EBF wake profile and certainly different 
than an axisymmetric jet profile. Since the wake induced detailed load 
distribution is determined by the wake velocity distribution, the use of 
a more realistic parabolic-type velocity profile inside the wake would 
affect the predicted flap load distribution. 

Finally, the wake boundaries are limited at the present time to 
elliptical shapes centered around the wake centerline. This boundary 
shape is reasonable for the forward part of the wake, but when substan- 
tial interaction with the wing or flap occurs, the boundary will distort 
and influence the widths and heights of the peak induced loadings. At 
the present time, the method cannot handle more complex shapes. 
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The predicted pressure distributions on the lifting surfaces in the 
region of maximum jet interference show the correct general trends, but 
the use of only a small number of discrete constant pressure panels 
prevents the prediction of the details of the pressures. Power-off 
pressure distributions are generally predicted quite well. 

The comparison of measured and predicted flow fields aft of EBF 
configurations show good similarity. Even though the two configurations 
are not the same, they have similar flap angles and identical thrust 
coefficients. The predicted flow field, corresponding to a converged 
solution, indicates similar flow angles and jet position beneath the 
wing as shown by the data. 

The preceding comments generally indicate that the overall momentum 
reaction and jet induced circulation effects are properly modeled, but 
the details of the loading distribution are not consistently adequate, 
based on the one set of data available. 

NIELSEN ENGINEERING & RESEARCH, INC. 

Mountain View, California 
November 1975 
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APPENDIX A 


INDUCED VELOCITY FIELD DUE TO AN 
ELLIPTIC CROSS-SECTION JET MODEL 


The concept of circular vortex rings is described in reference 14. 
If the vorticity is allowed to lie along an elliptical shape, the result 
is what shall be termed herein as an elliptic . vortex ring. The purpose 
of the following analysis is to determine the induced velocity field due 
to an elliptic vortex ring with constant vorticity and arbitrary 
eccentricity. 

Consider an ellipse specified by the equation 


Yi 

2 

a 


+ 



1 


(A-1) 


and a field point P located at the point as shown in fig- 

ure 6. A system of confocal elliptic coordinates (|,T]) is given by 


y = c cosh ^ cos t) 


z = c sinh ^ sin T) 


where 


= a 


- b-^ 


(A-2) 


(A-3) 


and is constant for all confocal ellipses. The semi-major and semi- 
minor axes, a and b, respectively, are given by 


a = c cosh 
b = c s inh 



(A-4) 


Therefore, from (A-2) 


Z = 
a 


cos 


b 


sin 


(A-5) 
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Note that constant values of ^ correspond to ellipses in which the 
distance between the foci is 2c, and constant values of rj correspond 
to hyperbolas confocal with one another and with the ellipses . Thus, 
and represent the ellipse and hyperbola 'through the point 

z^) and can be called the elliptic coordinates of P. 

Differential Induced Velocity Components 

The general formula for the induced velocity components at P(x^,y^, 
z^) due to a line segment of length ds and circulation r is given by 
the Biot-Savart law as 

dv = ("s X "r) (A-6) 

4ttR^ 


where 

R = vector from ds to P 

s = unit vector tangent to line segment; positive direction 
corresponds to direction of P 

Assxame the line segment ds to be in the plane x = 0. Therefore, 

"r = (Xq)"! + (y^ - y)'j + (z^ - z)lc (A-7) 


and 


s ds = dy j + dz k (A-8) 

where dy and dz lie on the ellipse described by equation (A-5) . 

Equation (A-8) can be written as 


s ds = ds = -a sin ri d-Q ] + b cos t) di] k 

Consider the confocal elliptic coordinates of P to be 

y = a cos n 
o o 'o 

z = b sin n 
o o ‘o 



(A-9) 


(-A-10) 
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The cross product in equation (A-6) can be carried out as follows, 
ds ("s X R) = d"^ X = 


X j k 

0 -a sin t] b cos ti 


y - y z - z 


= (ab - ab^sin rj^sin t] - a^b cos rj^cos T))dTi i 


+ (bz^cos T) df]) j + (az^sin t] dr|)k 


(A-11). 


From (A-6) we can now write the three components of the induced 
velocity in differential form. 


j, (ab - ab^sin 'H^sin t| - a^b cos rj^cos t]) dr^ 

dv = — — ■ 

X 4tt „3 


(A-12) 


p bx^cos T1 dT] 

■>''y - 4? ;3 


(A-13) 


_ ax sin ti dr] 

dv = 2 

z 4ir 3 

1\ 


(A-14) 


The radial distance R is obtained from (A-7) 


R^ = (x^) ^ + (a^cos - a cos t\)^ + (b^sin - b sin j])^ (A-15) 


Induced Velocity Components Due to First 
Quadrant of Elliptical vortex Ring 

Equations (A-12) , (A-13) , and (A-14) can be integrated over the 

first quadrant of the ellipse to obtain the portion of the induced 
velocities due to that part of the vortex ring. A variable change from 
1 ^ to t is made in accordance with the following transformation 
suggested in reference 18. 


t 


tan 


a 

2 


cos T] = 


1 - t 


1 + t" 


dr| 


2 

1 + t^ 


sin rj 


2t 

1 + t^ 


dt "1 




J 


(A-16) 
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Expanding in equation (A-15) results in 


= 


i “j t^ fc^cos^Ti + + b^ + + 2a a cos ti + 1 

(1 + t^)2 I L *o o o o 'o J 


+ t^[-4b b^sin ii^] + t^ [2 (c^cos^t)^ + b^ 

+ t[-4b b^sin 11^]+ [c^cos^t)^ + b^ + b^ + x^ - 2a a^cos + c^]} 

(A-17) 

Equation (A-17) can be rewritten as 

Tt + X«t + X t^ *4" X t + X 1 
L 3 2 1 oJ 


R^ = ^ 


where 


(1 + t^) 


Q = c^cos^Ti^ + b^ + b^ + x^ + 2a a cos n + c‘ 
'o o p o *0 


and 


X = 


X = 


X = 


c^cos^n + b^ + b^ + x^ - 2a a cos n + c 
1^0 o o o *0 

Q 

-4b b sin n_ 

o Jo 

Q 

2 (c^cos^tj^ + b^ + b^ + - c^) 


2 




X = X 
3 1 


J 


From equations (A-12) , (A-13) , and (A-14) , 


(A-18) 


(A-19) 


(A-20) 


dv = 

X 


(1 + t") 


2ttQ^/^ [t^ + X t^ + X t^ + X t + X 1^/^ 

+ t(-2ab sin Ti ) + (ab - a b cos n ) 
o 'o o 'o J 


[t^(ab + a^b cos t] ) 


(A- 21) 
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dv = 


2ttQ 


3/2 


bx^d - t^) (1 + t®) 


[f* + + X^t + X^J 


3/2 


dt 


(A-22) 


dv_ = 


ax t(l + t ) 
o 


27tQ 


3/2 


[t"* + Xgt^ + Xgt^ + X^t + 


(A-23) 


The limits of integration for the above velocity components are 
obtained from equation (A-5) . The lower limit at y = a, z = 0 
corresponds to r; = 0 which gives t = 0 from (A-16) . The upper limit 
at y = 0, z = b corresponds to r) = tt/ 2 and t = 1. Now define a 
set of integrals J^, to be 


J 


n 



t“^dt 


+ A. t + X t + 
2 1 



(A-24) 


The velocity components from (A-21) , (A-22) , and (A-23) become 


V 

X 


— 7 — f(ab + a b cos n ) J . - (2ab sin ti ) 

3/2 L o 'o' 4 ' o 'o 3 

2irQ 

+ 2ab J - (2ab sin ti )J, + (ab - a b cos ti ) J 1 
2 o 'o 1 o ‘o oj 


(A-25) 
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(A-26) 


(A-27) 


These velocity components are those induced by the first quadrant 
of the elliptic vortex ring. Details of the evaluation of the J 
integrals in equation (A-24) are presented in Appendix B. 


Induced Velocity Components Due to 
Complete Vortex Ring 

The velocity components induced by the four quadrants of an elliptic 
vortex ring are obtained from equations (A-25) , (A-26) , and (A-27) 

through an appropriate change in the coordinate system and a corresponding 
change in the sign of the resulting velocity component. Thus, the four 
components of the x velocity are 


5.7 


APPENDIX A 


X 
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= v„ at (x^.-yo^-^o^ 


X 


= V at (x^,y^,-z^) 


X 


(A-28) 


o’^o’ o' J 


where the subscripts I, II, III, and IV correspond to the particular 
quadrant of the vortex ring and is computed from equation (A-25) at 

the noted coordinates. 

The y velocity components are 


V. 


= V at (x^,y^,z^) 

V ' o o’ o 
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II 


= (^o’-^o’^o^ 
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III 
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IV 


= V at (Xo,-y^,-Zo) 


= Vy at (x„,yo»-^o) J 


(A-29) 


where v is computed from equation (A-26) at the noted coordinates. 

In the same fashion, the z velocity components are computed from 
equation (A-27) . 
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II 
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III 
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(A-30) 
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Each component of induced velocity at (x^,y^,z^) from a single ring 
is the sum of the contributions from the four quadrants. Some induced 
velocities from a family of elliptic rings with various axis ratios are 
presented in the results section of this report. 


Velocities Induced by Complete Jet Model 

The jet -wake model is made up of a series of vortex rings with their 
axes lying along the prescribed centerline of the jet as shown in fig- 
ure 7. If the jet is deflected spanwise under the wing, as shown in the 
top view of figure 7, the rings are assumed not to tilt by the centerline 
deflection angle, a. However, for vertical deflection, which is generally 
much larger, the rings tilt by the local angle 0, as shown in the lower 
view of figure 7. As the jet centerline is deflected, the orientation 
of a field point P with respect to a vortex ring changes with the ring 
position along the centerline. Since the jet induced velocities on the 
wing are required for interference calculations, the following set of 
equations transforms the induced velocities from the ring coordinate 
system to the wing system. 

The coordinates of P in the jet centered coordinate system are: 


"j. = 




> 


Z • = -Z., + Z- 

Dp P Q 


(A-31) 


where point Q denotes the center of the origin of the jet. Trans- 
forming (A-31) one step further, the coordinates of P in a ring 
centered coordinate system are 


x_ = (x . - X. ) cos 9 + (z. - z. ) sin 9 


= (y. - y^ ) 


z 


p 



X . ) sin 0 + (z . 

3s ^P 


z . ) cos 0 1 

3s J 


(A-32) 
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where the point S denotes the location of the center of a particular 
vortex ring. 

The three components of induced velocity in the ring system, v^, 

V . and V , are obtained by svimming the four parts of each component as 

y» z’ 

noted in equations (A-28) , (A-29) , and (A-30) , respectively. These 

velocity components, when resolved into components in the jet coordinate 
system, become 


u . = V cos 0 - V sin 9 
j X z 


V . = 

D 


V 

y 


> 


w. = V sin 0 + V cos 0 1 

■ j X z y 


(A-33) 


The induced velocities in the wing coordinate system are: 


u = -u . 
3 


V = V. V 


w = -w . J 


(A-34) 


The total velocity induced at P by the entire jet -wake model is 
obtained by summing the velocities induced by all the rings forming the 
wake . 
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EVALUATION OF J„ INTEGRALS FOR 
AN ELLIPTIC VORTEX RING 


The set of integrals 

j ^ f ^ 

^ J t* + X^t^ + X_t^ + X t + 




(B-1) 


appearing in equations (A-25) , (A-26) , and (A-27) need to he solved 

analytically in order to calculate the induced velocity components due 
to an elliptic vortex ring. 

The general solution of a quartic equation can be written as 
t"* + Xgt^ + X^t^ + X^t + X^ = (t - a) (t - a) (t - p) (t - 3) 


= [(t - b^)^ + a^][(t - b^)^ + aj] (b-2) 


where 


a = b 

+ ia , 

a = 

1 


S = b 

2 

+ ia , 

3 ’ 

3 = 


(B-3) 


The variable t in equation (B-2) is replaced by u using the 
transformation 


tn u = 


t - bi + a^g^ 


(B-4) 


from page 146 of reference 18. The constants in equation (B-4) are defined 
as follows. 


, a + a 


b = 

2 2 


1 4 


TTv 2 


a^ = 


2 _ 4a^ - (A - B)" 
(A + B)^ - 4a^ 


} 


(continued on next page) 


61 


APPENDIX B 


= (b - b + (a + a ) 
1 2 1 2 




= (b^ -b^)- + (a^ - a^)"-J 


From (B-4) 


t = 


(ai + g^b^)tn u + b^' - a^g^ 
1 + g^tn u 


which can be differentiated to give 


a (1 + g ) , 

. 1 ' dn u , 

dt du 

(1 + g^tn u) ^ cn^u 


Page 146 of reference 18 provides the relation 

dt 


[(t - a) (t - a) (t - p) (t - 3)] 


where 


"■/ JTT 


d(t> 


yj 1 - k^sin2 0 


= F(0,k) 


and 


g = 


A + B 


= 


4AB 


(A + B)' 


Siibstituting equation (B-6) into the first term of (B- 
in the following: 

2 2 (1 + g^) (1 + tn^u) 

(t - b ) + a = -i 

X 1 ,-1 . _ ..x2 


(1 + g^^tn u) 


(B-5) 

(concluded) 


(B-6) 


(B-7) 


(B-8) 


(B-9) 


(B-10) 


2) results 


(B-11) 
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Therefore, 

^ (1 + g^tn u)^ 

(t - a) (t - a) a^(l + g^) (1 + tn^u) 


(B-12) 


Applying the same svibstitution to the second term of equation (B-2) does 
not result in a simple relationship like equation (B-11) . The second 
term becomes 

(1 + g^tn u)^[(t - b^)^ + a^] = [ (a^ + b^g^)^ - 2b^g^ (a^ + g^b^) 

+ g^(a^ + b^)l tn^u 
12 2 J 

+ 2[(a^ + g^bj (b^ - a^g^) - b^ (a^ + b^g^) 
- b^g^(b^ - a^g^) + g^(aj + bj)] tn u 
+ - 2^2 

+ (a^ + b^)1 (B-13) 

2 S 

Ass\iming that (B-13) has the same form as (B-11) , the coefficient of the 
tn u term of (B-13) must be zero. This coefficient can be written as 

\ + [(^1 " ^2^^ + (^2 - ^i)]^! + (b^ “ ^ 2 ) = ° (B-14) 

which is a quadratic in g^ . Therefore, 

-[(b, -.?)]' +4a=(b3-b,) = 

^ 2a^ (b^ - b^) 

(B-15) 
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which can be rewritten in a general form and simplified as follows 


= 


-R + V + 


-R + V R^ + Q‘ 


-R + V 


R^ + 


J=Q- 


-R + 


V R^ + 


(B-16) 


Multiplying (B-15) by (B-16) results in an expression for g^ . 


[(b^ -bg)^ + (a| -ap] +^[(b^ -bj^ + (a^ + 4a^ ^ 

-[(^1 -^ J )] + Vf ^\ + 4 a ^(\ 


gf = 


(B-17) 


Substituting the relations for and B^ from equation (B-5) into 

(B-17) and. using the upper sign of the square root terms results in 

4a^ - (A - B)^ 

\ ^ i F 

^ (A + B)^ - 4a^ 


which is identical to the relationship in equation (A-5) . Thus, equa- 
tion (B-13) becomes 

(1 + g^tn u)^[(t - b^)^ + a^J = [^(a^ + b^g^)^ - 2 b 2 g^^ (a^ + b^g^) 

+ g^(a^ + b^)l tn^u 
12 2 J 

+ f(b, - a g ) ^ - 2b (b - a g ) 

L 1 1^1 21 11 

+ (a^ + b^)l (B-19) 

2 2 J 

A great deal of effort was expended to get the right-hand side of this 
equation into the same form as the numerator of equation (B-12) , but with 
no success. An alternate approach follows. 
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Substituting equation (B— 7) into (B-8) results in 


g = 


+ dn u 


(1 + g^tn u)^ cn^u ^ [^(t - b^)^ + a^J ^(t - b^)^ + a|J 


(B-20) 


or 


_ 2 / 1 1 \ 2 

a_ \ (1 + g^) 


[(t - + .J][(t - b,)» + aj] . (f) — 


_ dn^u 

4 4 

(1 + g tn u) on u 


(B-21) 


It can be shown that 


= (1 + tn^u) (1 + k'^tn^u) 

cn u 


(B-22) 


where 


k' = 1 - k" 


Equation (B-21) becomes 


af (1 + gf)^ 


f(t-b +an[(t-b +a^1 = i (1 + tn^u) (l+k'^tn^u) 

^ ^ ^ g^(l+g^tnu)^ 


(B-23) 


Removing the identity in (B-11) from this result produces the required 
relationship for the second term. 


[(i - b^)^ + aj] = 


(1 + gp (1 + k'^tn^u) 

g^ (1 + g^tn u) ® 


(B-24) 


Substituting equation (B-2) into (B-1) and letting n = 0, the 


integral for becomes 


1 

^ r dt 

i {[(t-b^)»][(t -b^,» 


(B-25) 
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Substituting equation (B-8) and (B-23) into (B-25) results in 




^ (1 + g^tn u)^ 

1 + ^ (1 + tn^u) (1 + k' 


^tn^u) 


du 


where the lower and upper limits are obtained from equation (B-4) 


(t = 0) 


(t = 1) 


tn 


tn 






1 - b^ + ^1?! 
\ 






The following relations are obtained from reference 18. 


1 + tn^u = ^ 


2 

cn u 




1 + k'Sn'u - (ffl-Jl) 

V cn u / ^ 

Siibstituting equation (B-28) into (B-26) gives 


:i + cs^)^ J 


a^d + g") 


/ T . \ cn u , 

(1 + g tn u) — — du 

1 j 2 

dn u 


A; 


g 

—2 / (cn u + g sn u) 

1 + J 


4 1 


a (1 + g ) 


dn^u 


du 


which, when expanded, results in the following integrals. 




• U 4 

A = I siL_u du 

O / ._2. 


dn u 


u. 


k ' k 


(1 + k'^)E(u) - 2k ' u - k (sn u) (cd u) 


U 




(B-26) 


(B-27) 


(B-28) 


(B-29) 


(B-30) 
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u. 


/ U 3 

du 

dn^u 


k"" L 


dn u + 


dn u 


-l^J 


u. 


■ U 2 2 

A = I dn 




L 


, 2 
dn u 


(1 + k'^)u - 2E(u) + k^ (sn u) (cd u) 




A 


u. 

-I 




k^ L 


U 3 

cn u sn u 
dn^u 


du 


-.u, 


dn u + 


dn u 


U 


A 


u. 

-J 


1 

k^ L 


U 4 
cn u 

, 2 

dn u 


du 


-2k*‘"u + (1 + k'^)E(u) - k^k'^(sn u) (cd u) 


A 


'"l 


where 



u = 

F(0, 

k)"" 

E(u) = 

E(0, 

k) 

\ 

sn 

u = 

sin 

/ 

0 

cn 

u = 

cos 

0 


(B-31) 


(B-32) 


(B-33) 


(B-34) 


(B-35) 
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Finally, 


2 v 2 


a^d + gp 


[A 4 + + 6 g^A^ + 4g^A^ + g^A^ (B-36) 


When n = 1 in equation (B-1) , 

1 


r = f — t dt 

^ J /r(t - b. )^ + a^l f(t - 


{[(t - hp^ + a^][(t - + a|]| 


21 ) 3/2 


(B-37) 


which can be simplified in the same manner described above for J^. 
Therefore, 


J i= 


a^(l + g^)^ 
1 • 


fU 

J (1 + g tn u)®[(a +b^g^)tn u + (b^ 


4 

v“l cn u ^ 


dn^u 


- a^g^) rU 


1 ~ ^ 3 1 

I (cn u + g, sn u) (cn u + t sn u) — ■ — du (B-38) 

a‘^(l + g 2)2 J ^ dn^u 

1 


where 


T = 




(B-39) 


Making use of the A^ integrals in equations (B-30) through (B-34) , 
equation (B-38) can be written as follows. 


•"l = .3 


g®(b^ - a^g^) 


a " (1 + g^)^ 

1 


[A 4 + (3g^ + t)A^ + 3g^ (g^ + t)A^ + (g^ + 3 t)A^ 


Similarly, 


(B-40) 


J = 


g^(b^ - a^g^)' 

a^d + g^)^ 

1 1 


[a^ + 2 (g^ + t)A^ + (gf + 4xg^ + t^) 


A 


+ 2xg^ (g^ + t)A^ + T^g^A^] 


(B-41) 
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g"^(b - a g 

J = [a + (3t + g, )A + 3T(g + t)A_ 

1 1 


+ T 


(t + 3g^)A^ + T^g^A ] 


g®(b - a,g,)^ . 2 3 4 -1 

j = [ A + 4 tA + 6t a + 4t a + t a 

4 2,t 2.2L4 3 2 1 o-J 

a (1 + g ) 

1 


(B-42) 


(B-43) 
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a/X 


0.500 

0.421 

0.354 

0.334 

0.323 

0.318 


0.525 

0.444 

0.376 

0,356 

0.345 

0.339 


0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

0.623 

0.959 

1.962 

3.600 

- 2.792 

- 1.263 

- 0.393 


0.535 

0.855 

1.829 

3.441 

- 2.936 

- 1.373 

- 0.460 


0.463 

0.770 

1.725 

3.321 

- 3.050 

- 1.467 

- 0.526 

- 0.222 

0.441 

0.746 

1.696 

3.375 

- 3.173 

- 1.495 

- 0.548 

- 0.240 

0.430 

0.733 

1.726 

5.571 

- 5.400 

- 1.509 

- 0.560 

- 0.250 

0.424 

0.728 

1.675 

3.265 

- 3.106 

- 1.516 

- 0.566 

- 0.255 


Table I.- Induced axial velocity in the plane 
of an elliptic vortex ring. 

















\z/b 

a/b\. 

0 

0.25 

0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

1.0 

0.358 

0.357 

0.346 

0.286 

0.204 

0.170 



0.004 

-0.035 

1.5 

0.326 

0.327 

0.321 

0.270 

0.192 

0.160 



-0.010 

-0.049 

3.0 

0.284 

0.286 

0.283 

0.238 

0.165 

0.133 

0.066 


-0.030 

-0.074 

5.0 

0.268 

0.270 

0.268 

0.224 

0.152 

0.120 

0.054 

0.023 

-0.048 

-0.086 

10.0 

0.259 

0.261 

0.259 

0.216 

0.144 

0.112 

0.046 

0.014 

-0.055 

-0.094 

00 

0.255 

0.257 

0.255 

0.211 

0.140 

0.108 

0.042 

0.010 

-0.060 

-0.098 



X 

b 


0.5, 


Y 

a 


0.0 


\z/b 

0 

0.25 

0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

a/b\^^ 











1.0 

0 

0.056 

0.129 

0.220 

0.259 

0.264 


0.242 

0.188 

0.102 

1.5 

0 

0.056 

0.130 

0.228 

0.276 

0.283 

0.278 

0.268 

0.214 

0.122 

3.0 

0 

0.055 

0.129 

0.231 

0.284 

0.293 

0.293 

0.284 

0.233 

0.140 

5.0 

0 

0.054 

0.128 

0.231 

0.235 

0.295 

0.296 

0.287 

0.238 

0.144 

10.0 

0 

0.054 

0.127 

0.230 

0.285 

0.295 

0.297 

0.288 

0.239 

0.146 

00 

0 

0.054 

0.127 

0.231 

0.286 

0.296 

0.297 

0.289 

0.240 

0.147 


vj 

■ to 


Table II.- Induced velocity field out of the plane 
of an elliptic vortex ring. 

















u^b 

r 


at 


X 

b 


0.0 


1 . 0 , 


Y. ^ 


a 


''V z/b 
a/b\ 

0 

0.25 

0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

1.0 

0.177 

0.170 

0.151 

0.118 

0.093 

0.085 

0.069 

0.061 


0.014 

1.5 

0.186 

0.180 

0.161 

0.127 

0.102 

0.093 

0.075 

0.066 


0.014 

3.0 

0.178 

0.172 

0.155 

0.123 

0.098 

0.089 

0.071 

0.062 

0.037 

0.007 

5.0 

0.169 

0.164 

0.147 

0.115 

0.091 

0.082 

0.064 

0.055 

0.031 

0.000 

10.0 

0.163 

0.158 

0.141 

0.109 

0.085 

0.076 

0.058 

0.049 

0.025 

-0.005 

00 

0.159 

0.154 

0.137 

0.106 

0.081 

0.073 

0.055 

0.046 

0.022 

-0.009 




w .b 

-f- at 

b 

»> i 

= 0.0 




N^^z/b 

a/n\ 

0 

0.25 

0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

1.0 

0 

0.033 

0.063 

0.084 

0.090 

0.091 

0.090 

0.089 

0.082 

0.064 

1.5 

0 

0.037 

0.071 

0.098 

0.106 

0.108 

0.108 

0.107 

0.100 

0.080 

3.0 

0 

0.039 

0.077 

0.108 

0.119 

0.121 

0.123 

0.123 

0.117 

0.098 

5.0 

0 

0.040 

0.078 

0.110 

0.122 

0.124 

0.126 

0.126 

0.121 

0.103 

10.0 

0 

0.040 

0.078 

0.110 

0.123 

0.125 

0.128 

0.128 

0.123 

0.105 

00 

0 

0.040 

0.078 

0.111 

0.123 

0.126 

0.128 

0.128 

0.124 

0.105 


Table II.- Continued 














































u^b 

r 


at 


X 

b 


1.5 


0.0 


Z = 


a 


\. z/b 
a/b\. 

0 

0.25 

0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

1.0 


0.082 

0.074 



0.049 


0.040 

0.031 

0.019 

1.5 


0.097 

0.088 



0.059 


0.048 

0.038 

:0.023 

3.0 


0.105 

0.095 

0.079 

0.068 

0.064 


0.053 

0.041 

0.025 

5.0 


0.102 

0.092 

0.076 

0.065 

0.062 

0.054 

0.050 

0.038 

0.022 

10.0 


0.098 

0.088 

0.072 

0.062 

0.058 

0.050 

0.046 

0.035 

0.018 

00 

0.098 

0 

0.085 

0.070 

0.059 

0.055 

0.047 

0.043 

0.032 

0.015 




w -b 

b 

i 

= 0.0 




\z/b 
a/b N. 

0 

0.25 

0.50 

0.75 

0.90 

0.95 

1.05 

1.10 

1.25 

1.50 

1.0 

0 

0.014 

0.027 





0.040 

0.039 

0.035 

1.5 

0 

0.018 

0.033 





0.051 

0.050 

0.046 

3.0 

0 

0.021 

0.040 

0.054 

0.060 

0.062 


0.064 

0.064 

0.060 

5.0 

0 

0.022 

0.042 

0.057 

0.063 

0.065 


0.068 

0.068 

0.065 

10.0 

0 

0.022 

0.042 

0.058 

0.064 

0.066 

0.068 

0.069 

0.070 

0.067 

00 

0 

0.022 

0.042 

0.058 

0.065 

0.066 

0.069 

0.070 

0.071 

0.067 


U1 


Table II.- Concluded 
































a 

(degs.) 


% 

0 

18.5 

— 

0.028 

2.3 

6.5 

0.14 



18.5 

1.10 

0.070 

4.0 

6.5 

0.54 

0.032 


18.5 

1.25 

0.074 


26.5 

1.72 

0.096 


Table III.- Combined lift due to leading-edge 
slat, fuselage, and nacelles; and pitching 
moment due to leading-edge slats on the 
EBF configuration of references 13 and 16. 


^T 

ram 

0.25 

0.057 

0.58 

0.088 

1.0 

0.118 


Table IV.- Ram drag coefficient for 
a single JT15D-1 turbofan 
engine, reference 13. 










(a) Wing coordinate systems. 

Figure 2.- Coordinate systems used in 
determining wing boundary condition. 




(a) Finite length vortex filament. 



(b) Semi-infinite vortex filament. 

Figure 3.- Influence function coordinate systems. 


80 


I 










I V 


Vj ,P(X,Y,Z) 



-*j AS I-*- 


Top view 
Side view 




Figure 7.- Vortex ring jet model. 


84 











Ref. 8 



Figure 9.- Measured half-width boundaries of circular 
elliptical, and rectangular cross-section jets 
exhausting into still air. 






Figure 10.- Spanwise variation of measured jet height 
distributions at approximately one flap chord aft 
of EBP model at various thrust coefficients. 



0 


2.0 3.0 4.0 


0.5 


1.0 

C 




(b) Equivalent radius ratio. 


Figure 11.- Measured jet expansion on 
an EBF model. 








a, degrees Cm 

(a) Lift and pitching moment coefficients. 


Figure 12.- Measured and predicted aerodynamic 
characteristics of an unswept wing with 
single slotted flap deflected 40°. 



wing x/c flap x/c 

(b) pressure distribution at t) = 0. 25. 

Figure 12.- Concluded. 
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Predicted Values 

Deflected 

Wake 

Undeflected 

Wcike 

0. 162 

0. 166 

-0. 168 

-0. 170 



Figure 13.- Measured and predicted pressure distribution 
at ri = 0.46 on a 35° swept wing with a 
partial span flap deflected 10°. 










(a) P/Pq = 1.0, y/a = 0.0 

Figure 16.- Induced axial velocity 
profiles inside elliptic jets. 








a 



0 0.2 0.4 0.6 0.8 1.0 1.2 

Uj/Y 


(c) P/Po = 1.0, Zj/b = 0.0. 
Figure 16.- Continued. 
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Eect 

:ted 

jet 


of jet expansion on 
axial velocity 
centerline. 

















Iteration 


Figure 21.- Convergence of the total wing-flap 
normal-force coefficient on a four-engine EBF 
configuration, = 0°/20°/40°, a= 18.5®. 



Iteration 


Figure 22.- Convergence of normal-force coefficients on the 

wing and flaps of a four-engir^ EBF configuration, 

6. = 0°/20°/40°, a = 18.5, C, = 4.0. 
r jx 
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0.2 


0.3 


0.4 


0.6 


0 


Figure 23.- 
coef f icients 


0. 5 

(a) Wing. 

Convergence of predicted spanwise distribution of section normal-force 
on a four-engine EBF configuration, 6f = 0°/20°/40°, C =4.0, a* 18.5°. 
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Initial jet —j 



Iteration 
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(b) Outboard jet, il = 0.42. 
Figure 24.- Concluded. 










I. 3 


(c) Flap 3. 


Figure 27.- Concluded 
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a, degrees 

(a) Lift and pitching-moment coefficients. 


Figure 28.- Measured and predicted longitudinal aerodynamic: characteristics 
of a four-engine EBF configuration, 6^ = 0°/20°/40°. 
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(b) Drag coefficient, 
igure 28.- Concluded. 



















(c) Flap 2. 


Figure 29.- Continued 









ref .16 




Theor 


Flap 3. 

- Concluded. 












(a) Wing. 


Figure 30.- Measured and predicted section normal force coefficients 


on the lifting surfaces of a four-engine EBF configuration 


6 ^ = 15°/35y55\ 


C . = 4.0 , a = 18.5°. 
M- 


I 
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0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


1 

(d) Flap 3. 


Figure 30.- Concluded. 
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(b) Flap 1. 


Figure 31.- Continued. 







Figure 31.- 
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T 
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Exp. (ref. 16) 
Theory 



edicted pressure coefficients on a four-engine 
5f = 15°/35°/55°/, C„ = 4.0°, a = 6.5°. 
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(a) Lift and pitching moment coefficient. 


Figure 33.- Measured and predicted longitudinal aerodynamic characteristics of 
a four-engine EBF configuration, 6^ =* 15°/35°/55°. 
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(b) Drag coefficient. 
Figure 33.- Concluded. 
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Figure 34.- Measured and predicted velocity profiles 
aft of EBF configurations, Tj « 0.25. 
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(a) 6^ = 17. 


Figure 35.- Measured and predicted longitudinal aerodynamic 
characteristics of a four-engine EBF configuration with 
an unswept wing and double- slotted trailing edge flaps. 
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Experiment, Ref. 17 
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Figure 35.- Concluded. 
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(a) Wing. 

Figure 36.- Effect of jet cross section shape on predicted spanwise distribution 
of section normal-force coefficients on the lifting surfaces of a four-engine 
EBF configuration, 5^ = 15°/35°/55°, C = 4.0, a = 18.5°. 
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(b) E 
Figure 36. 
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Figure 36.- C 
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